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Avant-Propos 



Ce memoire d'habilitation decrit I'application de procedures de renormalisation de 
type Ma-Dasgupta a plusieurs types de systemes desordonnes, classiques ou quantiques, 
dynamiques ou statiques. Afin de replacer dans une perspective commune ces differentes 
etudes, les premiers chapitres sont consacres a une presentation generale de la methode. 



Les chapitres generaux sur la methode 

• L 'introduction situe les renormalisations de type Ma-Dasgupta par rapport aux 
nombreuses approches qui existent pour les systemes desordonnes. 

Etant donne qu'une renormalisation de type Ma-Dasgupta est a la fois une "facon 
de penser" et une "fagon de calculer", ces deux aspects sont discutes separement, pour 
plus de clarte : 

• Le Chapitre n explique en details les idees physiques essentielles communes aux 
approches de type Ma-Dasgupta dans les differents contextes. 

• Le Chapitre [21 presente les regies quantitatives de renormalisation, qui donnent 
lieu a de jolis problemes de probabilites. 



Les chapitres sur des systemes desordonnes classiques 

La majorite du memoire est consacree a divers modeles de physique statistique 
classique, ou les approches de type Ma-Dasgupta sont utilisees soit pour caracteriser des 
proprietes de dynamique lente, de vieillissement ou de croissance de domaines, soit pour 
etudier les proprietes a I'equilibre thermodynamique de certains modeles desordonnes. 
Nous discuterons en particulier les modeles suivants : 

• la dynamique d'une particule dans un milieu aleatoire unidimensionnel, que ce 
soit dans un paysage de forces aleatoires (Chapitres |31 et 3)) , ou dans un paysage de 
pieges aleatoires (Chapitres El etEJ- 

• I'equilibre et la dynamique hors equilibre de croissance de domaines de chaines de 
spins desordonnees classiques (Chapitre 0), comme la chaine d'Ising en champ aleatoire 
ou la chaine verre de spin en champ magnetique exterieur. 

• la transition de delocalisation d'un polymere aleatoire, compose de monomeres 
hydrophobes et hydrophiles, a une interface entre deux solvants selectifs (Chapitre JHl). 
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Les chapitres sur des systemes desordonnes quantiques 

Apres ces modeles de physique statistique classique, nous reviendrons dans les deux 
derniers chapitres aux chaines de spins quantiques desordonnees, c'est a dire le 
domaine initial de la methode : 

• Le Chaoitre IIUI concerne la chaine de spin S = 1 antiferromagnetique aleatoire, 
qui presente une transition de phases en fonction du desordre, contrairement a la chaine 
S = 1/2. 

• Le Chapitre^^est consacre aux proprietes de taille finie de la transition de phase 
quantique de la chaine d'Ising avec couplages et champs transverses aleatoires. 

• Enfin, un petit apergu des differentes etudes de type Ma-Dasgupta qui existent 
dans la litterature, est donne en Annexe. 

Publications associees 

Ce memoire est base sur les publications notees [PI] ... [P15], publiees entre 1997 
et 2004, dont la liste se trouve a la fin, en page 11171 Par rapport au memoire, ces 
publications contiennent des discussions plus approfondies des modeles, des resultats 
complementaires, les derivations de tous les resultats, et des listes de references appro- 
priees. Par comparaison, le but de ce memoire est de donner une idee generale de la 
methode, des modeles qu'elle permet d'etudier, et des resultats qu'elle permet d'obtenir, 
sans entrer dans une description trop detaillee de chaque modele discute. 
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Introduction 



Ou Von situe les renormalisations de type Ma-Dasgupta parmi les approches 
qui cherchent a decrire les fluctuations spatiales dans les systemes desordonnes. 

La presence de desordre dans un systeme peut donner naissance a des phenomenes 
physiques tout a fait nouveaux, comme la localisation d'Anderson en matiere condensee, 
ou les proprietes de dynamique lente en physique statistique. L'etude des systemes 
desordonnes a done engendre un certain nombre d'approches specifiques depuis une 
cinquantaine d'annees : le premier exemple est la methode de mesure invariante de 
Dyson [HH] et Schmidt jl53j qui permet d'obtenir des resultats exacts pour les systemes 
unidimensionnels decrits par des produits infinis de matrices de transfert aleatoires 
J123lll5] . Pour comprendre I'interet et la specificite des procedures de renormalisation de 
type Ma-Dasgupta qui constituent le sujet de ce memoire, on peut classer les differentes 
facons d'aborder les systemes desordonnees en deux grandes categories : 

• il y a d'une part des approches qui commencent par moyenner sur le 
desordre, car elles se fixent pour objectif d'evaluer des observables automoyennantes, 
comme par exemple 1' energie libre s'il s'agit d'etudier la thermodynamique du systeme. 
II existe un certain nombre de prescriptions specifiques pour effectuer cette moyenne 
sur le desordre : la methode des repliques |133j . la methode supersymetrique ^\, la 
methode dynamique J23] (pour une presentation parallele des trois methodes : |117j V 
Apres cette moyenne sur le desordre, il n'y a plus d'heterogeneites spatiales, mais il y 
a en echange un systeme pur avec de nouvelles interactions effectives. 

• il y a d'autre part des approches qui cherchent a decrire les heterogeneites 
spatiales du desordre, comme certains arguments celebres et diverses procedures de 
renormalisation dans I'espace reel, que nous allons discuter un pen plus en details, car 
ces approches appartiennent a la meme 'famille de pensee' que les procedures de type 
Ma-Dasgupta. 



2 Introduction 

Des arguments probabilistes sur les fluctuations spatiales 
du desordre ... 

Par mi les arguments qui out joue un grand role dans la comprehension des systemes 
desordonnes, on pent citer 

(a) I'argument de Lifshitz J1191 I12UJ . qui permet de predire les singularites 
essentielles des densites d'etats en bord de spectre, en identifiant les configurations 
du desordre qui engendrent des etats dans cette zone d' energie, et en estimant les 
probabilites de presence de ces configurations favorables. 

(b) les phases de Griffiths , dans lesquelles des regions rares localement ordonnees 
induisent des singularites essentielles pour I'equilibre [HTj et des lois de decroissance 
lente pour la dynamique |148| 130] . 

(c) le critere de Harris |S9j sur la pertinence d'un faible desordre autour du 
point critique pur, qui estime I'influence des fluctuations locales du desordre sur les 
fluctuations de temperatures critiques. 

(d) I'argument de Imry-Ma jl()7j . qui permet de predire I'apparition de parois 
de domaines dans les systemes en champs aleatoires, en consider ant les fluctuations 
spatiales d' energie venant des champs aleatoires. 

(e) le theoreme de Chayes et al. |1J|, qui montre que les fluctuations spatiales 
du desordre impliquent I'inegalite u > 2/d pour I'exposant critique i' dans un systeme 
desordonne en dimension d. 

En fait, ces differents arguments ne font intervenir que deux proprietes statistiques 
difFerentes. D'une part, les arguments (a) et (b), qui sont tres proches |142( 1123] . re- 
posent tons les deux la prise en compte d'evenements rares : dans toute configuration 
infinie de desordre, il existe des domaines ordonnes arbitrairement grands avec des prob- 
abilites exponentiellement petites. D'autre part, les arguments (c), (d) et (e) reposent 
tons les trois sur le comportement typique en yiV de la somme d'un grand nombre A^ 
de variables aleatoires independantes. 

Ces arguments de type probabiliste sont bien fondes et presque "inattaquables" . A 
ma connaissance, le seul argument qui a donne lieu a une controverse jlU8j est I'ar- 
gument de Imry-Ma qui etait en desaccord avec la "reduction dimensionnelle" predite 
par les approches de theorie des champs, que ce soit en perturbation a tons les ordres 
Pl I1H4J ou dans le formalisme supersymetrique |146j . Les etudes rigoureuses |l()f)| IHH ] 
ont finalement donne raison a ... I'argument de Imry-Ma! Get exemple montre que 
les arguments dits "heuristiques" meme s'ils sont simples, peuvent avoir un contenu 
physique non-trivial, qu'il n'est pas toujours facile de retrouver par des methodes plus 
"sophistiquees" . 

En conclusion, ces arguments de type probabilistes permettent de bien comprendre 
la physique, car ils identifient les fluctuations locales du desordre qui sont responsables 
de tel ou tel phenomene. En contrepartie, il est souvent difficile de depasser I'aspect 
qualitatif et de 'calculer' vraiment... Pour aller au dela en restant dans le meme esprit, 
I'idee la plus naturelle est bien siir d'utiliser des procedures de renormalisation dans 
I'espace reel. 
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... aux procedures de renormalisation sur le desordre 

Le choix de travailler dans I'espace reel pour definir une procedure de renormalisa- 
tion, qui a deja un grand interet pour les systemes purs jl41llHl] . devient I'unique choix 
possible en presence de desordre si Ton souhaite decrire les heterogeneites spatiales. 

Renormalisations par blocs 

Les renormalisations par blocs, de type decimation ou de type Migdal-Kadanoff sont 
les renormalisations les plus utilisees pour les systemes desordonnes. Les procedures de 
type Migdal-Kadanoff J1351 llllj constituent en efFet des approximations simples pour 
effectuer des renormalisations par blocs sur les reseaux reguliers. Elles representent aussi 
des renormalisations exactes sur certains reseaux hierarchiques J112j . Parmi les systemes 
etudies, on pent citer par exemple le modele de Potts J116j . les ferromagnetiques dilues 
jlU9j . et surtout les verres de spin, qui ont donne lieu a un grand nombre de travaux : 
ces etudes concernent les diagrammes de phases |169( l27| l26] , diverses proprietes de la 
phase verre de spin |1621 11371 I152J . et surtout le caractere chaotique des trajectoires 
du flot de renormalisation J1311 1^1 11441 I162J qui est une grande nouveaute par rapport 
aux systemes purs. 

Par ailleurs, divers systemes desordonnes ont ete etudies par des procedures de 
renormalisation sur reseaux hierarchiques, en particulier le modele de Potts j5^ . les 
polymeres diriges en milieu aleatoire [SQiOESl, et les problemes de mouillage aleatoire 

inmni]- 

Citons enfin les procedures de renormalisation par blocs pour la diffusion dans les 
modeles de pieges jl24| et pour les chaines de spins quantiques aleatoires 92 , car nous 
discuterons dans ce memoire I'application des procedures de type Ma-Dasgupta a ces 
memes modeles. 

Renormalisation fonctionnelle pour les interfaces en milieux aleatoires 

Pour les modeles d'interfaces en milieux aleatoires, il existe une methode de renor- 
malisation fonctionnelle [ZHj : c'est une methode de theorie des champs qui etudie le 
flot de la fonction de correlation du desordre. Nous renvoyons a la revue |167j pour 
la description des divers developpements recents, et aux references I118J pour des 
comparaisons avec les methodes de repliques. 

Renormalisation pour les modeles XY desordonnes en 2D 

L 'introduction de desordre dans les modeles bidimensionnels de type XY, dans 
lesquels il existe des transitions de type Kosterlitz-Thouless dans le cas pur, a con- 
duit a une renormalisation de type gaz coulombiens, dans laquelle on etudie le flot 
d'une distribution de probabilite des fugacites jSJj afin de prendre en compte I'effet 
des heterogeneites spatiales sur les defauts topologiques. Cette approche permet aussi 
d'etudier la transition vitreuse d'une particule dans un potentiel aleatoire presentant 
des correlations logarithmiques [55] . 
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Renormalisations phenomenologiques pour les verres de spin 



La "theorie des droplets" 1741 pour les verres de spin en dimension finie, a pour orig- 
ine une renormalisation phenomenologique introduite par Mc Millan |129j et developpee 
par Bray et Moore [2H]- Dans la formulation de Bray et Moore ^^, I'idee essentielle est 
que la distribution de probabilite -Pl(J) des couplages effectifs J a I'echelle L converge 
vers une forme fixe, avec une largeur J{L) = JL^ qui depend de I'echelle L. L'exposant 
y et la distribution limite, qui sont calculables exactement en d = 1, ont ete etudies 
numeriquement en d = 2 (y < 0) et d = 3 (y > couplage fort ) [2H1- Cette idee d'un 
point fixe de couplage fort (ou de temperature nulle), decrit par une forme d'eclielle 
pour la distribution de probabilite d'une variable aleatoire, correspond en fait exacte- 
ment a la description que I'on obtient par les procedures de renormalisation de type 
Ma-Dasgupta lorsqu'on peut les utiliser. 

Renormalisation de Ma-Dasgupta pour les chaines de spins quantiques 

La procedure de renormalisation introduite par Ma-Dasgupta-Hu en 1979 |125j pour 
etudier la chaine de spin quantique 5 = 1/2 avec interactions antiferromagnetiques 
aleatoires a pour caracteristique essentielle de renormaliser de maniere inhomogene 
dans I'espace pour mieux s'adapter aux realisations locales du desordre. En effet, les 
methodes usuelles de renormalisation traitent I'espace de maniere homogene, en rem- 
plagant par exemple chaque bloc de spins de taille donnee par un super-spin. Si ce 
caractere homogene est naturel pour les systemes purs, on peut cependant se poser 
la question de sa legitimite pour les systemes desordonnes qui brisent I'invariance par 
translation du systeme. Ma, Dasgupta et Hu ont defini une procedure de renormalisa- 
tion sur r energie, et non sur la taille d'une cellule spatiale. La procedure consiste a 
decimer de maniere iterative les degres de liberte de plus haute energie, afin d'obtenir 
une theorie effective de basse energie pour la chaine de spin, qui porte aujourd'hui le 
nom de "random singlet phase". 

Cette procedure de renormalisation est en fait restee peu connue et peu utilisee 
pendant de nombreuses annees ... jusqu'aux travaux de Daniel Fisher en 1994-1995 
[751 [75] qui lui ont donne a la fois : 

(i) un statut theorique bien defini : les points fixes de "desordre infini" 

Alors que la methode apparaissait jusque la comme une procedure approximative 
peu controlee, Daniel Fisher a montre que le hot de renormalisation prenait une forme 
d'echelle qui convergeait vers un point fixe de "desordre infini" (ce qui signifie que le 
desordre est de plus en plus fort a grande echelle), ce qui rendait la methode asymp- 
totiquement exacte j7Sl|7n]. Par ailleurs, I'application de la methode a la chaine de 
spin avec couplages et champs transverses aleatoires ("Random Transverse Field Ising 
Chain" en anglais) lui a permis de montrer explicitement I'exactitude des resultats 
obtenus grace a une comparaison directe avec certaines observables calculees rigoureuse- 
ment pour le modele de McCoy et Wu |13L)1 1157] . (Ce modele desordonne de McCoy 
et Wu est un modele d'Ising classique bidimensionnel avec un desordre constant par 
colonne, qui est equivalent a la chaine RTFIC.) 
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(ii) des possibilites de calculs explicites assez remarquables 

Pour la chaine RTFIC [2^, la procedure de Ma-Dasgupta permet d'obtenir beau- 
coup de resultats nouveaux par rapport aux methodes rigoureuses |13Ul 11571 I143J . en 
particulier I'exposant critique exact /3 = (3 — V^)/2 pour raimantation spontanee. 
Plus surprenant, Daniel Fisher a meme calcule des observables que I'on ne connait pas 
pour le modele pur correspondant, c'est a dire pour le modele d'Ising pur bidimension- 
nel ! (par exemple la fonction d'echelle explicite decrivant I'ainiantation en fonction du 
champ applique dans la region critique) . 

Ces travaux de Daniel Fisher out engendre un grand interet pour ces methodes dans 
le domaine des spins quantiques (cf les Chapitres llOl et lTTl ainsi que 1' Annexe). Dans ce 
memoire, nous allons voir que les procedures de renormalisation de type Ma-Dasgupta 
ne se limitent pas aux systemes de spins quantiques, mais sont aussi un outil ideal pour 
etudier une classe beaucoup plus large de systemes desordonnes dans le domaine de la 
physique statistique. 
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Chapitre 1 

La signification physique des 
approches de type Ma-Dasgupta 



Oil Von explique la notion de "point fixe de fort desordre" 
et la nianiere de raisonner dans une approche de type Ma-Dasgupta 



Dans ce premier chapitre a caractere pedagogique, nous presentons en details les 
idees physiques essentielles communes aux diverses procedures de renormahsation utihsees 
dans divers contextes physiques. 



1.1 L'idee essentielle : le desordre qui domine par rapport 
aux fluctuations thermiques ou quantiques 

Dans tous les modeles consideres dans ce memoire, les differentes procedures de 
renormahsation reposent sur la meme idee : a grande echelle, le desordre domine par rap- 
port aux fluctuations thermiques ou quantiques. En particulier, les procedures de renor- 
mahsation de type Ma-Dasgupta sont vraiment specifiques aux systemes desordonnes 
et ne peuvent meme pas etre definies pour les systemes purs qui ne presentent pas 
d'heterogeites spatiales. 



1.1.1 Signification dans les differents contextes physiques 

D'une certaine maniere, les systemes purs, qui sont gouvernes par les fluctuations 
d'origine thermique ou quantique, sont caracterises par une grande "degenerescence" , 
dans la mesure oil tous les sites sont equivalents, alors que la presence de desordre 
leve completement cette degenerescence. Void quelques exemples de cette idee dans les 
differents domaines. 
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1.1.1.1 Exemple pour I'etat fondamental d'un systeme quantique 

Dans une chaine quantique antiferromagnetique pure de spin S = 1/2, on peut 
voir qualitativement I'etat fondamental comme une combinaison lineaire appropriee 
des etats qui correspondent a toutes les fagons possibles d'apparier les spins par deux 
pour former des singulets. En revanche, en presence de desordre, la procedure de renor- 
malisation de Ma-Dasgupta associe a chaque realisation du desordre un etat fondamen- 
tal qui correspond a une seule fagon d'apparier les spins en singulets. Ainsi, dans un 
echantillon desordonne fixe, un spin donne est completement correle a un seul autre 
spin de la chaine, qui peut eventuellement etre a une distance assez grande, mais n'est 
pratiquement pas correle avec les autres spins, meme ses voisins immediats sur la chaine. 

1.1.1.2 Exemple pour I'equilibre d'une chaine de spins classiques 

Dans la chaine d'Ising ferromagnetique pure, une parol de domaine peut se trouver 
de maniere equiprobable sur tons les liens car ils sont tons equivalents. L'argument 
classique energie/entropie entre 1' energie 2J que coiite une parol de domaine et I'en- 
tropie S ~ klnL associe a la position arbitraire de la parol dans un systeme fini de L 
sites, permet de comprendre I'absence d'ordre a longue portee a temperature finie et le 
comportement en Lt ~ e^'^' ^ de la taille typique d'un domaine. 

En revanche, la presence d'un champ aleatoire va lever cette equivalence entre 
tous les sites. L'argument d'Imry-Ma jlU7j . qui remplace l'argument energie/entropie 
precedent, est un argument energie/ energie : la comparaison entre 1' energie 2J que 
coiite une parol de domaine, et 1' energie typique ^/oL gagnee en profitant d'une fluc- 
tuation favorable de la somme X]j=i hi des champs aleatoires sur un domaine de taille 
L, conduit a une absence d'ordre a longue portee, meme a temperature nulle, et a une 
longueur typique Ljm ~ J'^/cr pour les domaines. Nous verrons que la procedure de 
renormalisation de type Imry-Ma permet de determiner les positions des parois des 
domaines d'Imry-Ma qui existent dans un echantillon donne. 

1.1.1.3 Exemples pour les marches aleatoires en milieux aleatoires 

Pour une marche aleatoire pure, la probabilite de presence se repartit de maniere 
gaussienne autour de I'origine, et 11 y a un etalement progressif au cours du temps. 

En revanche, pour une marche aleatoire dans un potentiel Brownien, 11 y a des 
regions particulieres a chaque echantillon qui localisent presque toute la probabilite 
de presence, et nous verrons comment les caracteriser par une procedure de type Ma- 
Dasgupta. En particulier, les fluctuations thermiques sont tout a fait sous-dominantes : 
la distance entre deux particules independantes ( c'est a dire deux histoires thermiques) 
qui diffusent dans le meme echantillon, reste une variable aleatoire finie dans la limite 
de temps infini, ce qui constitue la localisation de Golosov j85j . 

De meme, dans les modeles unidimensionnels de pieges caracterises par une distri- 
bution large des temps de piegeage p{t) ~ 1/t^+'^ avec < fi < 1, le front de diffusion 
dans un echantillon donne est essentiellement localise sur un nombre fini de pieges. A 
nouveau, nous verrons comment etudier leurs proprietes statistiques par une approche 
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de type Ma-Dasgupta generalisee. Au contraire, dans la phase ;U > 1 ou le temps moyen 
de piegeage est fini, il n'y a plus de localisation a temps infini, et une approche de type 
Ma-Dasgupta n'a plus de sens. 



1.1.1.4 Conclusion : les points fixes de fort desordre 

Si Ton s'interesse a la physique d'un systeme desordonne a grande echelle, il y a a 
priori trois possibilites pour revolution du desordre effectif par rapport aux fluctuations 
thermiques. En effet, lorsque I'echelle augmente, ce desordre effectif pent devenir 

(i) de plus en plus petit : le systeme est gouverne par un point fixe pur. 

(ii) de plus en plus grand : le systeme est gouverne par un point fixe de desordre 
infini. 

(iii) ou bien rester asymptotiquement a un niveau stable : le systeme est gouverne 
par un point fixe de desordre fini. 

Dans certains modeles, un desordre initial meme faible conduit a un point fixe 
de desordre infini (ii) : c'est en particulier le cas pour la chaine quantique antiferro- 
magnetique aleatoire S = 1/2, pour le modele de Sinai, pour la chaine d'Ising en champs 
aleatoires... 

Lorsqu'il existe un point fixe de desordre fini (iii), il est quantifie par un nombre 
qui varie continument, comme le parametre /U du modele de pieges. Le point fixe pent 
souvent etre qualifie de fort desordre dans une certaine region de parametres. Ainsi, 
dans le modele de pieges ou dans le modele de Sinai en presence d'un champ exterieur, 
la dynamique est gouvernee par un point fixe de fort desordre dans la phase /j, < 1 qui 
presente une localisation partielle du paquet thermique : il y a une probabilite finie 
que deux trajectoires thermiques dans le meme echantillon soient a une distance finie 
a temps infini. 

En conclusion, les methodes de renormalisation de type Ma-Dasgupta concernent : 

• les points fixes de desordre infini (ii). 

• les points fixes de desordre fini (iii) qui peuvent etre qualifies de desordre fort. 



1.1.2 Comment savoir si le desordre domine a grande echelle? 

1.1.2.1 Par des arguments theoriques a priori? 

L'importance relative qu'ont les fluctuations thermiques par rapport au desordre a 
grande echelle ne se voit pas directement sur le modele microscopique et n'est souvent 
pas tres bien connue pour la plupart des systemes desordonnes. Meme pour les systemes 
en champ aleatoire pour lesquels il existe un argument d'Imry-Ma energie/ energie 
discute ci-dessus, il n'existe pas, a ma connaissance, d'argument generalise qui inclurait 
dans la discussion les fluctuations thermiques des parois de domaines et qui estimerait 
l'importance de "I'entropie de decoupage en domaines Imry-Ma" dans un echantillon. 
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1.1.2.2 Par des etudes numeriques ? 

Dans les etudes numeriques qui ont a priori une information directe sur differentes 
configurations tfiermiques pour une realisation fixee du desordre, il est tres rare de 
trouver cette information sur un echantillon, car les resultats publies sont en general 
consacres a des moyennes de diverses observables sur les echantillons. Ce qui est bien siir 
dommage du point de vue des approches de type Ma-Dasgupta, ou I'information essen- 
tielle est justement I'importance des fluctuations thermiques dans un echantillon fixe. 
En effet, les moyennes sur les echantillons ont toujours le "risque" d'etre dominees par 
des evenements rares et de donner une fausse image des comportements typiques. Par 
exemple, dans le modele de Sinai, la largeur thermique moyennee sur les echantillons, 
qui est une observable naturelle a mesurer numeriquement pour caracteriser I'etalement 
du paquet thermique, diverge a grand temps. Ce resultat pourrait faire penser qu'il n'y 
a pas de localisation asymptotiquement, alors qu'en fait, la distance entre deux partic- 
ules independantes dans le meme echantillon reste une variable aleatoire finie a temps 
infini, ce qui correspond a une localisation tres forte en loi. 

1.1.2.3 Hypothese du desordre fort et sa verification 

En consequence, la demarche generalement utilisee dans les approches de type Ma- 
Dasgupta est la suivante : on commence par supposer que le desordre domine a grande 
echelle, et on verifie a la fin la consistance de I'hypothese. 

Plus precisement, les procedures de renormalisation de Ma-Dasgupta contiennent 
leur propre test de validite : si les distributions de probabilite ont une largeur qui croit 
indefiniment par le hot de renormalisation, elles conduisent a des resultats asympto- 
tiquement exacts, alors que si la largeur des distributions de probabilite converge vers 
une valeur finie, elles donnent des resultats qui seront seulement approches. 

1.2 La maniere de raisonner 

Dans la renormalisation de Ma-Dasgupta pour les chaines de spins quantiques, la 
procedure correspond a une decimation de degres de liberte sur I'Hamiltonien : elle 
est done assez proche des renormalisations usuelles, la seule difference etant que la 
decimation se fait de maniere iterative sur le couplage le plus fort, au lieu de se faire 
de maniere homogene sur toute la chaine a chaque etape. 

En revanche, dans les renormalisations de type Ma-Dasgupta pour les systemes 
de physique statistique qui seront discutes dans ce memoire, la maniere de raisonner 
s'ecarte beaucoup plus des procedures usuelles. En effet, le point de depart n'est pas 
une integration exacte ou approchee sur les degres de liberte du modele microscopique, 
c'est a dire sur la fonction de partition pour les problemes d'equilibre ou sur I'equation 
maitresse pour les problemes dynamiques. Le point de depart est plutot un argument 
physique heuristique, qui permet d'identifier a priori les degres de liberte qui sont im- 
portants a grande echelle. On definit alors directement la renormalisation sur ces degres 
de liberte juges importants, et on obtient ainsi dans les cas favorables des resultats ex- 



La maniere de raisonner 11 

acts dans la limite asymptotique ou la procedure de renormalisation est appliquee un 
grand nombre de fois. Comme ce melange entre arguments heuristiques au depart et 
resultats exacts a la fin peut paraitre deconcertant au premier abord, et se heurte meme 
souvent a une certaine incomprehension, il est utile d'analyser en details les differentes 
etapes du raisonnement sur le cas du modele de Sinai (ChapitreOl), qui est le cas le 
plus pedagogique. 

1.2.1 Modeles dynamiques : I'exemple du modele de Sinai 

La physique du modele de Sinai sera decrite dans le chapitre |31 qui lui est consacre. 
Notre but ici est uniquement de decrire la maniere de raisonner sur cet exemple precis. 

1.2.1.1 On identifie les degres de liberie du desordre qui vont etre impor- 
tants a grande echelle 

"On part d'un argument physique qualitatif : il existe depuis longtemps un argu- 
ment qualitatif simple H2] pour predire le comportement typique du deplacement en 
X r^ (Int)^ dans le modele de Sinai, au lieu du comportement habituel en x ~ \/t de 
la diffusion pure, que I'on peut resumer ainsi : le temps t{x) necessaire pour atteindre 
le point X > va etre domine par le facteur d'Arrhenius e^^ correspondant a la plus 
grande barriere B^ qu'il faudra franchir par activation thermique pour passer du point 
de depart x = au point x (Cette approximation par le facteur d'Arrhenius revient a 
effectuer une methode du col sur une expression analytique du temps du premier pas- 
sage). Dans un potentiel aleatoire Brownien, le comportement typique de la barriere 
est donne par Bx ~ ^/x, ce qui correspond a un temps d'Arrhenius t ~ e^^, ce qui 
correspond bien au scaling x ~ (Ini)^ apres inversion. 

"On prend au serieux cet argument pour determiner les degres de liberte qui vont 
etre importants a grande echelle" : I'argument heuristique ci-dessus suggere que les 
degres de liberte qui vont etre importants a grande echelle sont les grandes barrieres 
qui existent dans le potentiel aleatoire. Plus precisement, a un instant t fixe, la particule 
n'aura pas eu le temps de franchir par activation thermique les barrieres qui sont plus 
grandes qu'une echelle d'ordre (Tlni). 

1.2.1.2 On definit la renormalisation directement sur les degres de liberte 
importants du desordre 

On definit done une procedure de renormalisation sur les barrieres du potentiel 
aleatoire, dans laquelle on elimine de maniere iterative la plus petite barriere. Cette 
barriere minimum definit I'echelle de renormalisation T du paysage. Le paysage renor- 
malise a I'echelle T ne contient done que des barrieres plus grandes que F, toutes les 
barrieres plus petites ayant ete decimees. Lorsque I'echelle T augmente, la distribu- 
tion PriF) des barrieres F du paysage a I'echelle T tend vers une forme d'echelle 
9{F > r)P* ( p ), dans laquelle la distribution stationnaire P* caracterise le point 
fixe de desordre infini. 
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1.2.1.3 On etablit la correspondance avec le modele initial 

On associe a chaque temps t du modele initial, le paysage renormalise des barrieres 
a I'echelle F = Tint. On definit une dynamique effective sans fluctuations thermiques, 
dans laquelle la particule se trouve a I'instant t au minimum de la vallee renormalisee 
a I'echelle F = Tin t qui contient la position initiale a t = 0. On etablit aussi une corre- 
spondance entre les difFerentes observables physiques du modele initial et les proprietes, 
statiques ou dynamiques, du paysage renormalise. 



1.2.1.4 On verifie la consistance de la procedure dans le regime asympto- 
tique et on etudie les premieres corrections 

La probabilite que la particule ne soit pas dans la vallee renormalisee de la dy- 
namique effective, est d'ordre l/(lnt) et tend done vers zero dans la limite des temps 
infinis, ce qui montre la consistance de la procedure de renormalisation du paysage. 
On etudie les fluctuations thermiques par rapport a la dynamique effective, en con- 
siderant d'une part la distribution de probabilite a I'interieur de la vallee renormalisee, 
et d'autre part les evenements rares d'ordre l/(lnt) oil la particule n'est pas dans la 
vallee renormalisee de la dynamique effective. 



1.2.1.5 Discussion 

Get exemple montre bien comment cette facon de raisonner permet d'obtenir une 
description tres complete de la dynamique asymptotique a grand temps. II montre 
aussi tout I'interet de renormaliser de maniere inhomogene dans I'espace pour mieux 
s'adapter aux extrema locaux du desordre aux differentes echelles, par rapport aux 
renormalisations usuelles qui traitent I'espace de maniere homogene avec des cellules 
de taille fixee a chaque etape. 

D'une certaine maniere, la demarche que I'on vient de decrire utilise au maximum les 
idees qualitatives contenues dans la notion de renormalisation, comme la non-pertinence 
des details du modele microscopique sur les comportements a grande echelle, et la 
convergence vers des theories plus simples representant des classes d'universalite, avant 
de definir une procedure quantitative. II ne faut done pas reprocher aux approches 
de type Ma-Dasgupta de prendre comme point de depart des arguments physiques 
qualitatifs, car c'est de la justement que vient toute leur efficacite ! En effet, pour tous les 
modeles de physique statistique que nous allons considerer, alors qu'une renormalisation 
usuelle sur le modele microscopique desordonne n'aurait aucun espoir d'etre fermee et 
d'aboutir a la fin a des resultats exacts, I'approche de Ma-Dasgupta permet decrire une 
renormalisation fermee asymptotiquement exacte, directement sur les degres de liberte 
qui sont vraiment importants a grande echelle. 
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1.2.2 Modeles a I'equilibre : de I'argument Imry-Ma a la renormali- 
sation Ma-Dasgupta ! 

Pour les systemes desordonnes a I'equilibre qui seront discutes dans ce memoire, 
comme la chaine d'Ising classique en champ aleatoire (ou la chaine verre de spin en 
champ magnetique exterieur qui lui est equivalente par transformation de jauge), ou 
I'heteropolymere a une interface, la renormalisation de type Ma-Dasgupta est un moyen 
de construire les domaines Imry-Ma dans chaque echantillon. Les degres de liberte 
importants sont les sommes des energies aleatoires par domaines, et la renormalisation 
permet d'eliminer, de maniere iterative, les domaines qui sont instables vis a vis d'un 
retournement global, et d'obtenir ainsi, a la fin, la structure en domaines stables. 

1.3 Conclusion 

Le message essentiel de ce chapitre est done que les approches de type Ma-Dasgupta 
ont un sens chaque fois que les heterogeneites spatiales du desordre determinent I'etat 
dominant du systeme, alors que les fluctuations thermiques ou quantiques ne fournissent 
que des corrections sous-dominantes. Leur but est alors de construire I'etat dominant du 
systeme pour chaque realisation du desordre, generalement en fonction d'un parametre, 
qui est le temps pour les modeles dynamiques, ou la temperature pour les problemes 
d'equilibre thermodynamique. 

Pour mettre en oeuvre ce programme, il faut maintenant preciser les methodes de 
calculs utilisees pour effectuer quantitativement la renormalisation sur le desordre : 
c'est I'objet du Chapitre qui suit. 
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Chapitre 2 

Les regies de renormalisation de 
type Ma-Dasgupta 



Ou Von presente les regies quantitatives de renormalisation de type Ma-Dasgupta, 
en insistant sur la classe d'universalite des extrema Browniens. 



Alors que les renormalisations pour les systemes purs portent sur un nombre fini de 
constantes de couplages, les renormalisations pour les systemes desordonnes font inter- 
venir des distributions de probabilite, c'est a dire des fonctions qui appartiennent a un 
espace de dimension infinie, ce qui complique evidemment beaucoup I'analyse des flots 
de renormalisation et la recherche de points fixes. Cette difficulte conduit done le plus 
souvent a des etudes completement numeriques, ou alors, sur le plan analytique, a des 
approximations supplementaires qui consistent a projeter sur des espaces de dimension 
finie, c'est a dire a choisir une certaine forme analytique de distributions contenant 
quelques parametres dont on etudie les evolutions. Dans ce chapitre, nous allons voir 
que les regies de renormalisation de type Ma-Dasgupta engendrent, dans un certain 
nombre de cas favorables, des flots de renormalisation assez simples pour permettre 
un calcul exact de distributions de points fixes, qui ont souvent une interpretation 
probabiliste interessante. 

2.1 Qu'est-ce qu'une regie de renormalisation de type 
Ma-Dasgupta ? 

Si I'on souhaite englober les differentes renormalisations de type Ma-Dasgupta qui 
ont ete developpes dans les differents contextes physiques, on pent prendre comme 
definition generale la structure suivante : dans un systeme decrit par une collection 
de variables aleatoires {zi} qui ont des relations donnees de voisinage dans I'espace 
physique, la renormalisation consiste a decimer de maniere iterative la plus petite valeur 
min{zi}, avec eventuellement ses variables aleatoires voisines, afin de reconstruire une 
nouvelle collection de variables aleatoires {z^} reliees par de nouvelles relations de 
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voisinage dans I'espace physique. 

Au dela de cette definition generate un pen trop formelle, les deux caracteristiques 
essentielles sont les suivantes : 

• la renormalisation s'effectue sur la valeur extreme d'une variable aleatoire. Cette 
valeur extreme qui evolue par renormalisation constitue I'echelle de renormalisation : 
c'est le "cut-ofF' de la distribution renormalisee. 

• la renormalisation est locale dans I'espace : a chaque etape, il n'y a que le voisinage 
physique immediat de la variable aleatoire extreme qui est concerne par la renormali- 
sation. 

Nous allons maintenant decrire plus precisement les regies qui seront utiles dans les 
differents systemes desordonnes unidimensionnels discutes dans ce memoire. (Quelques 
exemples en dimension superieure seront evoques dans le Chapitre|Hlet dans I'Annexe.) 

2.2 La classe d'universalite des extrema Browniens 

Dans les systemes unidimensionnels avec desordre gele local, les approches de type 
Ma-Dasgupta conduisent souvent a une meme classe d'universalite qui a une interpretation 
probabiliste simple en termes d'extrema d'un potentiel Brownien a grande echelle. 

2.2.1 Les deux exemples quantiques fondamentaux 
2.2.1.1 La chaine antiferromagnetique aleatoire 



La regie de decimation proposee par Ma-Dasgupta-Hu |125j pour une chaine de spin 
S = 1/2 caracterisee par une suite { Jj} de couplages aleatoires positifs independants est 
la suivante : on cherche le couplage maximum dans toute la chaine, supposons pour 
fixer les notations que ce soit le couplage J2 = 0. On elimine ce couplage maximum 
et ses deux voisins, c'est a dire les trois couplages consecutifs {Ji,J2,J3), et on les 
remplace par un nouveau couplage efFectif 

f^^ (2.1) 

La signification quantique de cette regie sera discutee dans le Chapitre^l car ici notre 
but est de discuter uniquement I'aspect probabiliste de cette regie. Comme le nouveau 
couplage J' introduit est statistiquement independant de tons les autres couplages qui 
restent dans la chaine, cette regie de decimation engendre un flot de renormalisa- 
tion ferme pour la distribution Pn{J) des couplages J qui sont presents a I'echelle O 
de renormalisation, 011 O represente le couplage maximum qui evolue au cours de la 
renormalisation 



dn 



p{n, n) j dJa j dJb p{Ja, n) p{Jb, n)s(j- ^] (2.2) 
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2.2.1.2 La chaine d'Ising avec couplages et champs transverses aleatoires 
(RTFIC en anglais) 

Dans le modele RTFIC, il y une alternance le long de la chaine de couplages 
aleatoires Jj et de champs aleatoires hi, qui sont des variables aleatoires positives 
independantes, tirees a priori avec deux distributions differentes. La physique du modele 
sera discutee dans le chapitre ^2 Ici, notre but est simplement de decrire les regies de 
decimation [73] : on choisit a chaque etape le maximum des couplages et des champs 
Q = max{Ji, hj}. Si c'est le champ /i2 = fl, on I'elimine avec ses deux couplages voisins 
{J2, J3) pour former un nouveau couplage effectif 

J' = ^ (2.3) 

Si c'est le couplage J2 = ^, on I'elimine avec ses deux champs voisins (/ii,/i2) pour 
former un nouveau champ effectif 

, , All (lo , , 

h' = ^ (2.4) 

Ici encore, les nouveaux couplages introduits sont statistiquement independants de tons 
les autres couplages qui restent dans la chaine, ce qui permet d'ecrire un systeme ferme 
de deux Hots de renormalisation pour les distributions Pni-J) et Pc(/i) des couplages 
et champs effectifs a I'echelle Q. 

Nous allons maintenant expliquer comment ces regies qui apparaissent dans les 
chaines de spins quantiques peuvent s'interpreter en termes de statistique des extrema 
d'une marche aleatoire a grande echelle. 

2.2.2 Extrema des inarches aleatoires a grande echelle 

Dans plusieurs modeles de physique statistique que nous allons considerer, il existe 
un seul type de variables aleatoires (/,) sur une ligne, mais ces variables aleatoires 
peuvent etre positives ou negatives. Ces variables aleatoires representent par exemple 
les forces locales dans le modele de Sinai (Chapitres IHl et ^, les champs aleatoires 
dans la chaine d'Ising classique en champ aleatoire (Chapitre E}, les charges dans le 
modele du polymere hydrophile/hydrophobe (Chapitre inj. Dans les approches de type 
Ma-Dasgupta pour ces differents modeles, les degres de liberte importants du desordre 
sont les extrema du potentiel associe U{i) = YTj=o fj ^ grande echelle. 

Pour definir les extrema du potentiel initial, il faut done commencer par grouper 
ensemble toutes les forces fi consecutives de meme signe : on obtient alors un paysage 
constitue par une alternance de liens descendants {F^ , if) et de liens montants {F~ ,l~). 
Les barrieres F et les longueurs / sont maintenant des variables aleatoires positives 
qui representent respectivement les differences de potentiel et les distances entre deux 
extrema locaux consecutifs du modele initial. La regie de renormalisation du paysage 
est alors la suivante (Publications [P3,P4]) : on choisit la plus petite barriere T = 
min{Ff , F~}. Si la plus petite barriere est un lien descendant Ff = T, on I'elimine 
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F'= F,-F+F 

1 2 3 



ll I2 I3 



Fig. 2.1 - Illustration des regies de renormalisation pour un potentiel aleatoire unidi- 
mensionnel. 

avec ces deux liens montants voisins F^ et F2 pour former un nouveau lien montant 



(F-)' = Ff + F2 



(2.5) 



Si la plus petite barriere est un lien montant F-^ = T, on I'elimine avec ces deux liens 
descendants voisins F^ et F2 pour former un nouveau lien descendant 



{F 



+\' 



Fi+F, 



r 



(2.6) 



Ces regies sont done completement equivalentes aux regies de decimations (|2.3I2.4() du 
modele RTFIC par une simple transformation logarithmique sur les variables. Si Ton 
souhaite garder une information sur les distances initiales dans I'espace physique, il 
suffit d'ecrire la regie d'evolution de la longueur lors d'une decimation : la longueur du 
nouveau lien renormalise est simplement egale a la somme des trois liens elimines 



I' = li + l2 + h 



(2.7) 



Cette nouvelle longueur a encore la propriete d'etre independante statistiquement des 
autres longueurs qui restent dans la systeme. En revanche, la longueur va etre correlee 
avec la barriere F qui existe sur le meme lien. Ainsi, on peur ecrire un systme ferme 
de deux equations de flots pour les deux lois jointes P^{F, I). 



2.2.3 Distribution des extrema du potentiel Brownien pur 

Le cas symetrique fi = des modeles statistiques correspond au point critique 
quantique du modele RTFIC In J = In/i ou I'aimantation spontanee s'annule. A grande 
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echelle, les liens descendants et les liens montants deviennent statistiquement equivalents, 
et le seul parametre physique pertinent est la variance ff = 2a du modele initial, comme 
dans le theoreme de la Limite Centrale. II existe alors un point fixe de "desordre infini" 
caracterise par une distribution jointe P*(ry, A) des variables d'echelle rj = -^— pour 
les barrieres et A = ^^ pour les longueurs. Cette loi jointe qui caracterise le point fixe, 
est definie par sa transformee de Laplace par rapport a la longueur A [73] 

[ °^ dXe-P^V*{r],X) = 9{r,>0) ^ ^-vVp^o^^Vp (2.8) 

Jo smh ^ 

En particulier, la distribution des barrieres seules est une simple exponentielle 

V*i7]) = 0iv>0)e'~'^ (2.9) 

alors que la distribution des longueurs seules a la forme d'une serie infinie d'exponen- 
tielles 

-^EMr(-i)-'"^'''* 



(les deux series se correspondent par une formule d'inversion de Poisson). 

La convergence vers la solution de point fixe H2.8() est d'ordre l/L pour les marches 
aleatoires |75) . En revanche, pour le mouvement Brownien unidimensionnel qui represente 
deja la limite continue universelle des marches aleatoires, le point fixe ()2.8|) est un 
resultat exact a toute echelle F, comme le montre un calcul direct par des integrales de 
chemin contraintes (Publication [P9]). 

2.2.4 Distribution des extrema du potentiel Brownien biaise 

Dans le cas dissymetrique ori il existe un biais fi = fo> 0, les equations de renor- 
malisation conduisent a une famille de solutions a un parametre note 5 |75j : les deux 
distributions jointes Pp {F, I) des liens descendants et montants out pour transformees 
de Laplace 



r^ dle-^^P*^{F,l) = e{F > r) A/^^^^!g^,-(^-r)[V^coth v^^^.] ^^.u) 
JO sinh ^/p + S^ 

En particulier, les distributions de barrieres seules ont les formes exponentielles suiv- 
antes 

n-i^) =''(^>r)3?^»-'"'"^,^-„Jt»-'^-^'=^ (212) 

P*-{F) =e{F>r) _^l2^r e~^^~^^T3JW ^ OiF > r)26e-^^-^^^^i2.13) 
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Du point de vue du paysage renormalise, la signification du parametre 26 est done 
claire : la distribution P^~ (F) des grandes barrieres opposees au biais /o reste stable 
a grande echelle (en dehors de la presence du cut-off F) et (25) est le coefficient de 
la decroissance asymptotique exponentielle de cette distribution. II faut maintenant 
preciser la signification du parametre 6 par rapport au modele microscopique de depart. 
Pres du point critique correspondant k 6 = 0, le parametre 5 pent etre developpe 
au premier ordre dans le biais [TS] 

5 = - + 0(/o) (2.14) 

a 

Si Ton s'ecarte du voisinage immediat du point critique, la bonne definition non- 
perturbative du parametre 6 en terme de la distribution initiale Q{f) des variables 
(fi) du modele microscopique de depart est que 6 represente la solution de I'equation 
(Publication [P4]) 

r+oo 

e-25/, = / ^fQ^fy-2Sf ^ -^ (2.15) 



Cette definition correspond bien siir exactement, au changement de notation pres 25 = 
HJT ^ a la definition du parametre sans dimension /i qui represente I'exposant de la phase 
de diffusion anormale x ^ t^ des modeles de Sinai avec biais J1141 IS71 1^ , dont nous 
reparlerons dans le chapitre 0] Le developpement dans les deux premiers cumulants 
redonne bien I'expression simple (|2.14jl . qui n'est exacte a tons les ordres que pour une 
distribution gaussienne. 

A nouveau, si Ton considere directement le mouvement Brownien unidimensionnel 
biaise comme limite continue des marches aleatoires biaisees, les solutions ()2.11() sont 
des resultats exacts a toute echelle F, comme le montre un calcul direct par des integrales 
de chemin contraintes (Publication [P9]). 

Alors que la solution (|2.8j) du cas symetrique 5 = est un point fixe de "desordre 
infini" , la solution H2.13() du cas dissymetrique 5 > est un point fixe de "desordre fini" , 
car la distribution Pp^ (F) des grandes barrieres opposees au biais a asymptotiquement 
une largeur finie ^ . La renormalisation usuelle de type Ma-Dasgupta ne donne done des 
resultats exacts que dans la limite 5 — > 0. Nous discuterons en details dans le Chapitre 
|1] comment generaliser la methode usuelle lorsque le parametre b est petit. 

2.3 Renormalisation d'un potentiel aleatoire unidimen- 
sionnel quelconque 

Comme la renormalisation de type Ma-Dasgupta consiste a etudier les extrema d'un 
potentiel unidimensionnel lorsqu'on ne garde que les barrieres superieures a une certaine 
echelle F, il est possible de definir la procedure pour un paysage aleatoire quelconque. 
En particulier, la procedure de renormalisation pent etre implementee numeriquement 
pour des potentiels correles, ce qui a ete fait pour le cas de correlations logarithmiques 
[59j . La Publication [P9] etudie ce que Ton pent dire de maniere analytique pour le cas 
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des potentiels Markoviens. Nous decrirons en particulier la solution explicite pour le 
cas d'un potentiel Brownien en presence d'un potentiel deterministe harmonique dans 
le Chapitrel?! 

2.4 Regies de decimation plus generales 

Dans certains systemes unidimensionnels, en particulier pour la chaine quantique 
de spin S = 1 discutee dans le Chapitre ^1 on est amene a definir des procedures 
de renormalisation plus compliquees que la renormalisation d'un potentiel unidimen- 
sionnel, ce qui peut donner lieu a des phenomenes plus riches : en particulier, pour la 
chaine quantique de spin S = 1, cela permet de decrire une transition de phase de type 
percolation a temperature nulle. 

2.5 Variables auxiliaires et exposants critiques 

Dans les renormalisations de type Ma-Dasgupta, on appelle 'variables auxiliaires' les 
variables associees aux liens qui vont evoluer selon des regies de decimation paralleles 
lorsque leurs variables principales associees sont renormalisees. Le premier exemple 
important est la longueur I que nous avons deja rencontree (|2.7j) . C'est une variable 
auxiliaire car ce n'est pas elle qui determine la renormalisation, mais la barriere associee 
F : en effet a chaque etape, on ne choisit pas la longueur minimale, on choisit la barriere 
F minimale. 

Une autre variable auxiliaire importante pour le RTFIC est I'aimantation m des 
amas de spins : cette variable n'existe qu'en association avec les champs aleatoires hi 
et evolue en association avec la regie ()2.4() selon [7S] 

m' = mi + 7T1-2 (2-16) 

Plus generalement, dans les differents modeles discutes dans ce memoire, que ce soit 
pour la diffusion de Sinai (Chapitre OJ, pour la chaine d'Ising classique en champ 
aleatoire (Chapitre EJ, et pour les problemes de reaction-diffusion (Publication [P5]), 
on est conduit a etudier diverses variables auxiliaires qui evoluent selon la regie generate 

m' = anil + bm2 + cms (2-17) 

oil (a, b, c) sont des constantes. 

Nous avons deja vu comment I'etude jointe des barrieres et des longueurs car- 
acterisent la statique du paysage renormalise a une echelle donnee H2.8|) : en particulier, 
dans le paysage renormalise a I'echelle P, dans lequel il ne reste que des barrieres F > T, 
les longueurs ont pour scaling I ~ P^, ce qui est le scaling Brownien usuel comme il 
se doit. En revanche, en dehors de ce cas tres particulier a = 6 = c = 1, les variables 
auxiliaires ()2.17|) conduisent a des exposants non-triviaux m ~ P* qui reffetent des pro- 
prietes 'dynamiques' de la renormalisation sur toutes les echelles precedentes P' < P, et 
qui contiennent plus d' informations que le paysage a I'echelle P seulement. Par exemple 
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pour le RTFIC, la regie tres simple (|2.16|) correspondant aa = c = let6 = conduit 
pour raimantation ?n, ~ F''^ a I'exposant irrationnel egal au nombre d'or [^ 

^(a = c = l,6 = 0) = i±^ (2.18) 

Get exemple montre qu'une simple marche aleatoire unidimensionnelle contient des ex- 
posants non-triviaux si Ton s'interesse a des proprietes un peu 'subtiles' qui concernent 
la dynamique du paysage des extrema en fonction de I'echelle. 

Plus generalement, pour une variable auxiliaire de type (|2.17j) . on obtient que, 
lorsque la condition a + c = 2 est satisfaite (ce qui arrive assez souvent en pratique pour 
les observables physiques les plus naturelles), I'exposant associe satisfait une equation 
du second degre et vaut 



0(<. + c = 2,6) = l±4±5 (2.19) 

En revanche, lorsque a + c 7^ 2, I'exposant (f){a, b, c) satisfait une equation plus com- 
pliquee faisant intervenir la fonction hypergeometrique confluente U{A, B, z) (Publica- 
tions [P4] et [P5]). 

2.6 Comparaison avec certains modeles de croissance 

Les decimations de type Ma-Dasgupta ont une grande parente avec certains modeles 
de croissance, qui ont ete introduits de maniere completement independante. Ces modeles 
geometriques de croissance considerent une suite d'intervalles sur une ligne qui evolue 
par une transformation iterative sur le plus petit segment qui reste, avec les diverses 
regies suivantes : 

(i) dans le "cut-in-two model" j61j . I'intervalle le plus petit est elimine et donne 
une moitie de sa longueur a chacun de ces deux voisins. Ce modele introduit done des 
correlations entre les intervalles voisins et a ete etudie numeriquement j61j . 

(ii) dans le "paste-all model" [OJ, I'intervalle le plus petit est elimine et donne 
toute sa longueur a I'un de ces deux voisins tire au hasard de maniere equiprobable. 
Ce modele n'introduit pas de correlations entre les intervalles voisins, et la distribution 
invariante des longueurs a ete calculee [HJ . 

(iii) dans le "instantaneous collapse model" |151j . I'intervalle le plus petit est elimine 
avec ses deux voisins pour former un seul nouveau domaine I' = li + 12 + h- Ce modele 
decrit en fait la dynamique effective a grand temps d'un champs scalaire unidimension- 
nel qui evolue selon une equation de Ginzburg-Landau a temperature nulle J14UI I151J , et 
a done suscite un grand interet en tant que modele soluble de croissance. Les resultats 
exacts concernent la distribution invariante des longueurs jlSlj ). I'exposant de persis- 
tance |31) qui caracterise la variable auxiliaire d' = di+d^, I'exposant d'autocorrelation 
[SU qui caracterise la variable auxiliaire q' = qi — q2 + 93 et enfin I'exposant de persis- 
tance generalise |126) qui caracterise la variable auxiliaire m' = mi + pm2 + m^ avec 
un parametre p. 
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La seule difference 'technique' avec les regies de decimation de type Ma-Dasgupta 
est que, dans les modeles de croissance, c'est la longueur qui est la variable principale 
qui determine la renormalisation, alors que dans les modeles desordonnes, la longueur 
n'est qu'une variable auxiliaire, la variable principale qui definit la dynamique etant 
une variable de desordre. Cette difference explique les differences analytiques entre les 
solutions de points fixes et d'exposants dans les deux types de modeles. 

D'un point de vue physique, cet exemple montre qu'une dynamique sans desordre 
intrinseque, definie a partir d'une condition initiale aleatoire, pent etre gouvernee, d'une 
certaine maniere, par un 'point fixe de desordre infini'. 

2.7 Conclusion 

La specificite des regies de renormalisation de type Ma-Dasgupta est de decimer 
localement une variable extreme de desordre de maniere iterative. Cette structure tres 
particuliere permet souvent d'obtenir en dimension d = 1 des solutions explicites pour 
les distributions de probabilite du paysage et pour les exposants critiques associees aux 
variables auxiliaires representant des observables physiques interessantes. 

Ce Chapitre 2 termine la presentation generale des methodes de type Ma-Dasgupta. 
Le reste du memoire est consacre a I'application de ces methodes a divers modeles 
desordonnes de physique statistique ou de matiere condensee. 
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Chapitre 3 

Marche aleatoire dans un 
potentiel Brownien 

3.1 Presentation du modele 

Le modele d'une marche aleatoire dans un potentiel aleatoire Brownien, qui porte 
aujourd'hui le nom de "modele de Sinai", a beaucoup interesse les mathematiciens 
probabilistes depuis les travaux de Solomon jl58j Kesten et al. |114j et Sinai jl56j , ainsi 
que les physiciens des systemes desordonnes depuis les travaux initiaux de Alexander 
et a/.m et Derrida-Pomeau [SZl- Depuis, il y a eu de nombreux developpements dans 
les deux communautes : la revue recente |171j contient une presentation des references 
importantes pour les mathematiciens probabilistes, alors que les differentes approches 
des physiciens sont presentees dans les revues (^01 1^ I2H • 

Pour les physiciens, I'interet du modele de Sinai est double : d'une part, la marche 
de Sinai represente un modele dynamique simple en presence de desordre, sur lequel 
on pent tester un certain nombre de concepts generaux, et d'autre part, elle apparait 
naturellement dans divers contextes, par exemple dans la dynamique d'une parol de do- 
maine dans la chaine d'Ising en champs aleatoires (cf ChapitreEl) ou dans la dynamique 
d'ouverture de la double helice d'ADN en presence d'une force exterieure |122j . 

3.1.1 Version continue : Diffusion dans un potentiel Brownien 

La version continue du modele de Sinai correspond a I'equation de Langevin j22j 

- = -U'{x{t)) + v{t) (3.1) 

dans laquelle r]{t) represente le bruit thermique usuel 

< v{t)v{t') >= '^T5{t - t') (3.2) 

et U{x) est un potentiel aleatoire Brownien 



{U{x)-U{y)Y = 2a\x-y\ (3.3) 
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Plus generalement, dans tout le memoire, la moyenne thermique d'une observable / 
sera notee < / >, et la moyenne sur le desordre d'une observable / sera notee /. 

3.1.2 Version discrete : Marche aleatoire sur reseau ID 

Dans la version discrete sur reseau f Figure ITTj) . la particule qui se trouve sur le 
site i a une probabilite cui de sauter vers la droite et une probabilite {1 — coi) de sauter 
vers la gauche. Les Wj sont des variables aleatoires independantes dans ]0, 1[. La marche 
aleatoire est recurrente seulement si Inwj = ln(l — uji), ce qui constitue le cas de Sinai, 
et qui correspond a I'absence de biais du potentiel U{x) (|3.3|) de la version continue. 



1-co. 




i-l i i+l 

Fig. 3.1 - Version discrete du modele de Sinai 

3.2 Principe de I'approche de renormalisation 

Nous avons deja explique dans le premier chapitre (Page [TT|) la maniere de raisonner 
sur I'exemple du modele de Sinai. 

L'idee essentielle de notre approche est de decomposer le processus xu,r){t)-, representant 
la position de la marche aleatoire engendree par le bruit thermique r/(t) dans le potentiel 
aleatoire Brownien U{x), selon une somme de deux termes 

X{u,v} (t) = iTT'iu} (t) + y{UM (t) (3.4) 

• Le processus rn^ijj{t) porte le nom de "dynamique effective" dans les 
publications et represente la position la plus probable de la particule a I'instant t : elle 
correspond au meilleur minimum local du potentiel aleatoire U{x) que la particule a 
eu typiquement le temps d'atteindre depuis sa condition initiale pendant I'intervalle de 
temps t. Comme le franchissement d'une barriere de potentiel F necessite un temps 
d'Arrhenius d'ordre tp = tqc^^ , on pent etudier en detail cette dynamique effective en 
utilisant une procedure de renormalisation de type Ma-Dasgupta, qui consiste a decimer 
de maniere iterative les plus petites barrieres qui existent dans le systeme. On associe 
alors au temps t le paysage renormalise dans lequel il ne reste plus que des barrieres 
plus grandes que I'echelle de renormalisation F = Tint, ce qui correspond a une echelle 
F^ = (Tlnt)^ pour les longueurs. La position mfuj{t) correspond alors au minimum 
de la vallee renormalisee a I'echelle F = Tint qui contient le point initial. 

• Le processus y{u,r)}{'t) represente I'ecart par rapport a la dynamique effective. 
Dans la limite de temps infini, c'est un variable aleatoire qui reste finie, ce qui con- 
stitue le phenomene de localisation de Golosov : toutes les particules qui diffusent 



Etude de la dynamique effective 
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Fig. 3.2 ~ Figure tiree de la reference 140^ de J. Chave and E. Guitter, avec sa legende 
"Evolution with time (in logarithmic scale) of the distribution P{x, t) in a given energy 
landscape (drawn below). The evolution runs over lO''' iterations. The intinsity in the 
grey scale is proportional to (— In P{x, t)), i.e. darker regions correspond to higher values 
ofP{x,t)." 

dans le meme echantillon a partir du meme point de depart avec des bruits ther- 
miques r] differents sont asyniptotiquement concentrees dans la meme vallee de mini- 
mum ?n,|f/}(t). Plus precisement, si on considere les premieres corrections a grand temps, 
la probabilite qu'une particule ne soit pas dans la vallee correspondant a la dynamique 
effective m^ijy{t) est d'ordre l/(lnt), auquel cas la particule se trouve a une distance 
d'ordre (Int)^ de mu{t). Ces evenements sont done rares (leur probabilite tend vers 
zero a grand temps) mais ils dominent cependant certaines observables, comme par 
exemple la largeur thermique Ax'^{t) ~ < y^(i) > ~ (Int)^ qui diverge. 



3.3 Etude de la dynamique effective 



3.3.1 Front de diffusion 



A cause du phenomene de localisation de Golosov, la distribution de la variable 
d'echelle X = tttttttt, par rapport au bruit ry dans un echantillon fixe est asympto- 

-^ — est la variable 



(Tlnt)2 

tiquement une distribution delta de Dirac 6{X 



M) ou M 



(Tlnt)^ 

d'echelle de la dynamique effective. Si Ton souhaite maintenant calculer la moyenne du 
front de diffusion moyen sur I'ensemble des echantillons (ou I'ensemble des conditions 
initiales, ce qui revient ici au meme), il suffit d'etudier la distribution de M, et on 
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F'= F,-F+F 

1 2 3 



ll I2 I3 



Fig. 3.3 - Principe de la renormalisation pour etudier la dynamique effective m{t) 



retrouve ainsi tres simplement la loi de Kesten 



P{X) = LT 



-1 



p^\X\ 



p \ cosh ^ 



-TT ^-^ 



vr ^ 2n + 1 

n=0 



'\\n 2 

iy_g-^{2n+l)|X| 



(3.5) 



qui est un resultat exact des mathematiciens probabilistes J115| IHH]. Get example montre 
explicitement comment la procedure de renormalisation de type Ma-Dasgupta permet 
d'obtenir des resultats asymptotiques exacts, et donne confiance dans les resultats nou- 
veaux qu'elle donne pour des proprietes plus fines pour lesquelles il n'existe pas de 
"theoreme" correspondant. 



3.3.2 Distribution de 1' energie 

De meme, la variable d'echelle pour 1' energie 

[/(x(0)) - U{x{t)) U{m{0)) - U{m{t)) 



w 



(Tint) 



(Tint) 



est entierement determinee a grand temps par la seule dynamique effective. Cette vari- 
able reduite a pour loi limite quand t — > 00 : 

V{w) = e(w <l){A-2w- 4e""') + e{w > 1) (2e - 4) e""' 

La loi est continue, ainsi que sa derivee en w = 1, mais la derivee seconde est dis- 
continue en w = 1, ce qui peut sembler surprenant! En effet, pour tout temps fini, la 
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D(w) 




w 



Fig. 3.4 - Distribution limite de la variable d'echelle pour l' energie. 



distribution de 1' energie est analytique, et ce n'est que dans la limite de temps infini 
que la loi limite de la variable d'echelle presente une discontinuite de la derivee seconde. 
II est interessant de noter que dans le travail recent d'un mathematicien probabiliste 
[51] sur le temps de retour a I'origine apres I'instant t, il apparait aussi une distribu- 
tion asymptotique non- analytique pour une variable d'echelle qui represente aussi une 
variable de type energie et qui correspond aussi au point w = 1 dans nos notations. 



Nous avons aussi calcule la loi limite jointe de la position X 



energie w 



U{x{0))-U{x(t)) 



x(t)-x{0) 
(Tlnt)2 



et de 



(Tint) 

transformation de Laplace 



qui est definie par les deux expressions suivantes en 



+00 



-.-sX 



dXe-'^V{X,w > 1) 



sinh-y/s 



gv^coth Vi _ 2 cosh ^\ e-^'v^coth V3 



+00 



dXe 



-sX 



ViX w <l) = s™'^\^(^~"^) _ sinh2^ ^_^^/^^„th^/i 



(3.6) 



3.3.3 Proprietes de vieillissement 



Le front de diffusion a deux temps P(x, t;x', t'|0,0) presente un regime de vieil- 
lissement en (Int/lnt'). Dans les variables d'echelle X = (x/ln t) et X' = (x'/ln t), 
le front de diffusion est determine par la dynamique effective. La procedure de renor- 
malisation permet de calculer la loi jointe des positions {m(i),m,(t^)} a deux temps 
successifs t >tw (Publication [P4] ). En particulier, ce front de diffusion a deux temps 
presente une fonction 5{X — X') de Dirac, qui traduit le fait qu'une particule peut 
se trouver piegee dans une vallee dont elle n'arrive pas a sortir entre t' et t. Le poids 
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D{t, t^) de cette fonction S{X — X') qui represente la probabilite d'avoir m{t) = m{tw) 
a pour expression 

D{t,t^) = I (^^y (5-2e^-(^)) (3.7) 

Ce resultat caracterise bien ce qu'on appelle "vieillissement" : plus tw est grand, plus 
la particule a trouve une "bonne" vallee et plus il faudra de temps pour en sortir. 



3.3.4 Statistique des retours a I'origine de m{t) 

La procedure de renormalisation permet aussi de montrer que la distribution de la 
sequence Fi = Tlnti, r2 = Tlnt2 ••• des temps de retour a I'origine de la dynamique 
effective m(t) a une structure simple : c'est un processus Markovien multiplicatif defini 
par la recurrence r^+i = OkTk, dans laquelle les coefficients {oj} sont des variables 
aleatoires independantes, de loi 

, , 1 / 1 1 \ , 3±^/5 , , 

^(") = 7! (^iTA- - ^TTatJ ^-^^ ^± = -^- (3.8) 

Deux consequences import antes : 

• Le nombre total R{t) de retours a I'origine pendant [0, t] se comporte en 

i?(t)~ -ln(rini) (3.9) 

o 

alors que le nombre total S{t) de sauts pendant [0, t] se comporte en 

5(t) ~ -ln(rint) (3.10) 

o 

(mais ici il y a des correlations entre les temps de sauts.) 

• La probabilite que m{T) > pour r g]0, t] met en jeu un exposant de persistance 
irrationnel 



n(t) 



1 



(Tlnt)2 



— 3 — \/^ 
avec 6 = ^— = 0.19... 



alors que la probabilite qu'un marcheur donne x{t) ne repasse pas par I'origine x(o) 
pendant ]0,t] a un exposant de persistance simple 



liiit) 



(Tlnt)2 



avec 9 = — 
2 
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3.4 Etude des proprietes de localisation 

3.4.1 Distribution asymptotique du paquet thermique 

La distribution de la position relative y = x{t) — m{t) par rapport a la position m[t) 
de la dynamique effective a pour limite a temps infini une distribution de Boltzmann 
dans une vallee Brownienne infinie 

ou Ton moyenne sur deux trajectoires Browniennes {Ui,U2} qui forment une vallee 
infinie. Cette formulation est bien equivalente au theoreme de Golosov [SSj. La loi 
P{y) pent etre calculee en transformee de Laplace en termes de fonctions de Bessel 
(Publication [P8]) En particulier, on obtient le comportement asymptotique algebrique 

Piy) ~ 472 (3-12) 

qui pent etre interprete ainsi : alors que les configurations C/(y) typiques donnent une 
decroissance en e~^v^ pour le facteur de Boltzmann, il existe des configurations rares 
qui reviennent pres de f7 ~ a une grande distance y avec une probabilite en l/(y^'^). 
Nous avons aussi calcule la fonction de correlation a deux particules 

/ + 00 

dx{P{x, t|xo, 0)P(x + /, t|xo, 0)] (3.13) 

-oo 

qui decroit aussi algebriquement en 1//^". 

3.4.2 Parametres de localisation 

Les parametres de localisation, qui mesurent les probabilites moyennes de trouver 
k particules au meme point a temps infini 

n= lim / dx[P(x,t|xo,0)]'= = — ^(a/32)'=-i (3.14) 

*^°°J-oo r(2A;) 

ont un comportement a grand k qui est domine par les vallees tres etroites qui ont une 
petite fonction de partition (Publication [P8]). 

3.4.3 Comparaison avec les fonctions d'equilibre 

Ces differentes observables qui caracterisent la statistique du paquet thermique dans 
la diffusion de Sinai coincident en fait avec leurs analogues statiques definies par la limite 
thermodynamique de la distribution de Boltzmann dans un potentiel Brownien sur un 
intervalle (Publication [P8]). Cette convergence vers I'equilibre du paquet thermique 
(alors que la dynamique effective reste indefiniment hors equilibre) n'est plus vraie des 
qu'on ajoute une force constante (cf Chapitre|^ 



32 



Marche aleatoire dans un potentiel Brownien 



3.4.4 Largeur thermique et evenements rares 







(d) 



Fig. 3.5 - Representation des trois types (a) (h) (c) d'evenements rares d'ordre 1/T 
qui dominent la largeur du paquet thermique a grand temps ( la croix represente la 
condition intiale). Le cas (d) est un exemple d'evenement rare d'ordre l/F^ qui donne 
une contribution sous-dominante. 

La decroissance algebrique (|3.12|) de la distribution asymptotique de la position 
relative implique que le second moment < y^ > diverge a temps infini. Pour obtenir son 
comportement a grand temps, il faut prendre en compte les evenements rares suivants 
(Publication [P6]) : 

(a) une vallee renormalisee pent avoir deux minima eloignes dans I'espace qui sont 
presque degeneres en energie. 

(b) deux barrieres voisines peuvent etre presque degenerees. 

(c) une barriere pent etre en train d'etre decimee (F + e) 

Ces evenements rares, representes sur la Figure 13.51 out tons les trois une proba- 
bilite faible d'ordre 1/F, mais ils donnent lieu a division du paquet thermique en deux 
sous-paquets, separes par une grande distance d'ordre F^. En consequence, ce sont ces 
evenements qui dominent la largeur thermique 



< x^{t) > 



<x(t)>2 oc -(F 

t^oo i 



2n2 



r(rinT)= 



(3.15) 



et plus generalement tons les moments divergents d'ordre k > 1/2 



\x{t)- < x{t) > 



■ ~ cfeT(rinr) 



2fc-l 



(3.16) 



Les constantes c^ en prefacteur peuvent etre calculees a partir des proprietes statistiques 
des evenements rares (a,b,c) decrits ci-dessus (Publication [P4]). 
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3.5 Relation avec la theorie generale des dynamiques 
lentes et des etats metastables 

3.5.1 Etats metastables 

Dans la description qualitative usuelle des dynamiques lentes, que ce soit pour les 
verres, les milieux granulaires ou les systemes desordonnes, la notion d'etat metastable 
joue un grand role. Comme les etats metastables au sens strict n'existent que dans les 
approximations de champ moyen ou dans la limite de temperature nulle, si on veut 
utiliser ce concept pour des systemes en dimension finie a temperature finie, il faut 
considerer des etats metastables de temps de vie finie [20], en separant la dynamique 
en deux echelles de temps : il y a d'une part des degres de liberte "rapides", qui 
atteignent vite un quasi-equilibre local, ce qui correspond aux "etats metastables", et 
il y a d'autre part une dynamique lente hors equilibre, qui correspond a 1' evolution des 
etats metastables. 

Dans ce langage, notre description de la marche aleatoire de Sinai pent etre refor- 
mulee ainsi : 

• Les etats metastables a I'instant t sont les vallees du paysage renormalise a Techelle 
r = Tint : en effet, les marcheurs qui sont partis a t = d'un point de cette vallee 
n'ont pas eu le temps d'en sortir a I'instant t. 

• Dans chaque vallee renormalisee, il y a un quasi-equilibre decrit par une distribu- 
tion de Boltzmann a I'interieur de la vallee. 

• La dynamique lente correspond a revolution du paysage renormalise avec Techelle 
r = Tin t : certains etats metastables disparaissent et sont absorbes par un voisin. 

3.5.2 Observables a un temps et Conjecture d'Edwards 

Comme la 'Conjecture d'Edwards', qui propose de calculer les quantites dynamiques 
par une moyenne plate sur les etats metastables, a donne lieu a beaucoup de travaux 
recents ^2 ESI) il est interessant de reconsiderer de ce point de vue notre procedure. 

Dans notre approche, toutes les observables a un temps se calculent effectivement 
par une moyenne sur les vallees renormalisees (qui sont les etats metastables), mais 
avec une mesure qui depend de 1' observable : 

• Pour une condition initiale uniforme, la taille du bassin d'attraction d'une vallee 
est donnee par sa longueur, et done on utilise une mesure ponderee par la longueur 
J dllP{l) pour calculer la distribution de la position, de 1' energie, etc... 

• Pour les observables du paquet thermique, nous avons utilise une mesure plate 
sur les vallees Browniennes infinies, car il y a une independance entre la taille de la 
vallee et la statistique du fond de la vallee. 

3.5.3 Decomposition du front de diffusion sur les etats metastables 

Si Ton veut decrire a la fois la dynamique effective des vallees et I'equilibre de 
Boltzmann dans chaque vallee renormalisee, on pent ecrire le front de diffusion dans 
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un echantillon sous la forme 

P{xt\xoO) c^ V ^e-^^(^)ey,(x)%(:Eo) (3.17) 

La sonime porte sur toutes les vallees renormalisees Vp a I'echelle F = T In t. La notation 
9v{x) designe la fonction caracteristique de la vallee V, c'est a dire 6v{x) = 1 si a; 
appartient a la vallee et 9v{x) = sinon. Enfin Zy = Jydxe~^^^' represente la 
fonction de partition de la vallee V . 

Cette expression du front de diffusion (|3.17() correspond tout a fait a la construction 
generale en presence d'etats metastables (cf 20l I16UJ et les references incluses), dans 
laquelle I'operateur d'evolution e~^^^ est remplace par un projecteur sur les etats (i) 
d'energie Ei < 1/t 

e-*^^^~^|P, xQil (3.18) 

i 

L 'interpretation est claire : "everything fast has happened and everything slow has not 
taken place" |16Uj . Pour le modele de Sinai, les expressions explicites sont les suivantes : 
les etats a droite 

sont bien positifs, normalises, avec des supports qui ne se recouvrent pas, alors que les 
etats a gauche 

Qi{x) = e{x e V^'^) (3.20) 

sont simplement egaux a 1 sur le support de leur vecteur propre a droite associe, et 
nuls ailleurs. 

Dans le modele de Sinai, on pent en fait aller au dela de cette description a un 
temps en terme de projection sur les etats metastables, en considerant la dynamique 
de ces etats metastables pour obtenir des informations sur les proprietes spectrales de 
I'operateur de Fokker-Planck. 



3.6 Etude des fonctions propres de I'operateur de Fokker- 
Planck 

3.6.1 Construction des fonctions propres de I'operateur de Fokker- 
Planck associe a la dynamique effective 

II est interessant de voir comment change I'expression (|3.17() lors de la decimation 
d'une vallee renormalisee, c'est a dire lors de la disparition d'un etat metastable. La 
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decomposition de I'operateur d'evolution sur les valeurs propres £^„ > et les fonctions 
propres a droite <I>^ et a gauche #^ de I'operateur de Fokker-Planck Hpp s'ecrit 

P{xt\xoO) =< x\e~^"pp \xo >= Y^ e-^"^^^{x)^^{xo) (3.21) 

n 

En dehors de I'etat fondamental n = d'energie nuhe £"0 = qui correspond a 
I'equihbre de Boltzmann sur I'echantillon total 



$^(x) = l/^/z^t (3.22) 

<(x) = e-^(^)/^/yZto; (3.23) 

la comparaison avec I'equation (|3.17|) conduit aux identifications suivantes pour les 
etats excites n > 1 : les energies En sont determinees par les echelles de renormalisation 
r„ = Tlnt„ = — Tlnii^„ qui correspondent aux decimations de barrieres. Lors d'une 
decimation, deux vallees Vi et V2 se joignent en une seule nouvelle vallee renormalisee 
V', et les fonctions propres associees s'ecrivent 



^^(^) = J /\^P (^gvi(^) - -^Ov,{x)) (3.24) 

$«(x) = e-^(^)/^$^(a;) (3.25) 

On pent verifier que ces fonctions propres ont toutes les proprietes requises d'orthonor- 
malisation 

j dx^^{x)^^{x) = 5n,m (3.26) 

et de normalisation de la probabilite totale du paquet thermique 

dx^^{x) = (3.27) 



Au dela du modele de Sinai, la structure 1)3. 25(1 en termes de fonctions de partition 
partielles semble done decrire plus generalement les fonctions propres de I'operateur 
de Fokker-Planck pour les dynamiques lentes dans lesquelles les etats metastables dis- 
paraissent de maniere hierarchique et emboitee. 

3.6.2 Structure spatiale des fonctions propres : 2 pics et 3 echelles 
de longueur ! 

Un etat propre (|3.25|) presente deux pics qui correspondent aux minima des vallees 
Vi et V2. Chacun des deux pics a une largeur finie, qui represente la longueur car- 
acteristique associe poids de Boltzmann a temperature T autour du minimum d'une 
vallee. La distance entre les deux pics est d'ordre 1{E) ~ F^ ~ (Ini?)^. Loin des min- 
ima, mais a I'interieur de la vallee renormalisee r < F^, la fonction propre ^n{x) a 
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une decroissance gouvernee par le poids de Bolzmann e~^^^'^> de comportement typ- 
ique e~'^^. En particulier au bord de la vallee r ~ F^, Tainplitude typique est d'or- 
dre e~^^. Au dela des deux vallees en jeu, Tapproximation simple 1)3. 25() avec des 
fonctions theta devient insuffisante. Pour estimer la decroissance d'un etat propre sur 
une distance r > F^, il faut considerer II39J que les deux points sont separes par un 
nombre de vallees renormalisees d'ordre ^7 st que le recouvrement entre deux vallees 
voisines n'est pas zero, mais d'ordre e~'^ ^. Une theorie de perturbation conduit alors a 
une decroissance asymptotique exponentielle en e"'^ ^ , ce qui correspond en effet a la 
longueur de localisation X{E) r^ T r^ (—T\nE) calculee exactement par une methode 
de type Dyson- Schmidt pour le probleme de Schrodinger associe [22 . 

En conclusion, les proprietes des fonctions propres font done intervenir trois echelles 
de longueur qui coexistent : 

• une echelle finie I ~ 1 qui caracterise la largeur d'un pic, et qui est reliee a la 
localisation de Golosov du paquet thermique. 

• une echelle 1{E) ~ (Inii^)^ qui represente la distance entre les deux pics et qui est 
reliee a la distance totale parcourue au temps t ~ 1/E. 

• une echelle X{E) ~ (— InE) qui caracterise la decroissance exponentielle asympto- 
tique de la fonction d'onde, et qui correspond a la longueur de localisation du probleme 
de Schrodinger associe. 

3.7 Conclusion 

Le modele de Sinai est un exemple parfait de 'point fixe de desordre infini'. La 
procedure de renormalisation donne une image tres complete de la dynamique asymp- 
totique. Elle permet d'obtenir des resultats exacts explicites sur la dynamique effective 
d'une particule, sur les proprietes de vieillissement, sur les proprietes de localisation 
du paquet thermique, et sur les evenements rares qui gouvernent la largeur thermique. 
Par ailleurs, la description du modele de Sinai en termes de vallees renormalisees est 
un exemple explicite de la theorie generate des etats metastables de temps de vie fini, 
et permet d'obtenir une image tres claire de la structure des fonctions propres de 
I'operateur de Fokker-Planck. 

Mentionnons pour terminer que certains resultats obtenus par notre approche de 
renormalisation de type Ma-Dasgupta ont ete depuis confirmes par des etudes de 
mathematiciens probabilistes, que ce soit sur le poids de la partie singuliere du front de 
diffusion a deux temps [SS] ou sur la statistique des retours a I'origine de la dynamique 
effective ^. 

Publications associees 

• Principe de la renormalisation et premiers resultats : Publication [P3] 

• Dynamique effective, evenements rares et vieillissement : Publication [P4] 

• Sur la dynamique de I'energie : Publication [PIO] 

• Sur la localisation de Golosov et I'operateur de Fokker-Planck : Publication [P8] 



Chapitre 4 

Marche aleatoire dans un 
potentiel Brownien biaise 

4.1 Presentation du modele 

4.1.1 La phase de diffusion anormale 

L'introduction d'une force constante Fo dans I'equation de Langevin (|3.H) du modele 
de Sinai discute au Chapitre precedent est evidemment tres naturelle. Ce modele avec 
force a meme encore plus interesse les mathematiciens et physiciens depuis longtemps, 
car il presente une serie de transitions de phase dynamiques |114[ 1571 \'2'Z\ en fonction du 
parametre sans dimension fj, = FoT/a. En particulier, il existe une phase de diffusion 
anormale < /i < 1 qui est caracterisee par le comportement asymptotique sous-lineaire 



< xit) > ~ t'' (4.1) 



alors que pour // > 1, la vitesse devient finie : < x{t) > ~ V{fi)t. 

4.1.2 Le modele dirige de pieges 

II a ete propose depuis longtemps un argument heuristique j^ H2] selon lequel le 
modele de Sinai biaise etait asymptotiquement equivalent a un modele dirige de pieges 
defini par 1' equation maitresse 

dPtin) PAn) PAn - I) , , 

— --^ = i^^ + ^^^ L (4.2) 

di T„ 7:„_i 

dans laquelle les r„ sont des variables aleatoires independantes distribuees avec une loi 
algebrique 

9(^) ^ -TfH (4-3) 

La phase de diffusion anormale < /i < 1 correspond alors a la phase dans laquelle la 
moyenne des temps de piegeage est infini. 
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4.2 Principe de la renormalisation 

La renormalisation presentee dans le ChapitreElpour la marche de Sinai pent etre 
appliquee en presence d'un biais, mais elle ne donne des resultats exacts dans la limite de 
temps infini t — > oo que si le biais tend vers zero /i — > : par exemple, la renormalisation 
donne comme front de diffusion une loi exponentielle pour la variable d'echelle X = ^^ 
ce qui ne coincide avec le resultat exact mettant en jeu une loi de Levy J1141 [^ que 
dans la limite ^ — > 0. Le fait que la dynamique effective ne soit plus exacte lorsque n est 
fini a pour origine le fait que la distribution des barrieres -F_ opposees au biais H2.13() ne 
devient pas une distribution de largeur infinie par renormalisation, mais qu'elle converge 
vers une distribution exponentielle de largeur finie proportionnelle a 1/(2(5) = T/n. 
Ceci montre que la propriete de concentration du paquet thermique tout entier dans 
la meme vallee renormalisee a grand temps, qui est valide dans la limite ;U ^ 0, et 
qui constitue la localisation de Golosov, n'est plus exacte pour // fini. Nous avons done 
propose de generaliser la methode de renormalisation en incluant un certain etalement 
du paquet thermique sur plusieurs vallees renormalisees. Nous allons d'abord discuter 
cette procedure generalisee pour le modele dirige de pieges : dans un echantillon donne, 
la distribution de probabilite d'un marcheur est une somme de distributions 5 de Dirac, 
avec la structure hierarchique decrite dans la legende de la figure (|4.H) . 
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Fig. 4.1 - Structure des pieges importants pour une particule partie de x = a t = : 
La ligne pointillee separe les "petits" pieges (qui ont un temps de piegeage Ti < t) et 
les "grand" pieges (qui ont un temps de piegeage Ti > t). Le premier grand piege note 
M est le piege principal, qui est occupe avec un poids d'ordre 0(1). Le grand piege 
suivant Li et le plus grand piege Si parmi les petits pieges situes avant M sont les 
pieges secondaires occupes avec un poids d'ordre 0{fi). Le troisieme grand piege L2, le 
plus grand piege I2 parmi les petits pieges situes entre M et Li, et le second plus grand 
piege S2 parmi les petits pieges situes avant M sont les pieges tertiaires occupes avec 
un poids d'ordre 0(^^). 
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4.3 Resultats pour le modele dirige de pieges 




10000 
n 



15000 



20000 



Fig. 4.2 - Reproduction d'une figure de la reference \4^ de A. Compte et J. P. Bouchaud 
avec sa description : "Distribution of probability after a time t = 7x 10"'^'^ for a particular 
sample of disorder in our ID directed random walk model with /i = 0.4. The simulation 
was done with 1.000 particles in a lattice o/ 20.000 sites. This probability distribution is 
made of several sharp peaks that gather a finite fraction of the particles. However, the 
position of these peaks is scattered on a region of space of width f^. As time progresses, 
the position and relative weights of these peaks of course change, but at any given (large) 
time only a finite number of peaks, corresponding to very large trapping times, contain 
most of the particles. " 

A partir de cette description du front de diffusion dans chaque ecliantillon, il est 
possible de calculer exactement des series perturbatives en fj. pour toutes les observables. 
En particulier, les calculs explicites jusqu'a I'ordre ^^ (Publication [Pll]) du front de 
diffusion de la variable d'eclielle X = -^j de la largeur thermique 



lim 

i"»oo 



< An2(t) > 



IT- 



f2/^ 



^(2 In 2) + //2[ + 2 In 2(ln 2 - 2 + 2je)] + 0(/i 

6 



et du parametre de localisation 

Y2{fi) = l-^(21n2)+/x2(41n2 



vr 
~6 



+ 0{fi' 



(4.4) 



(4.5) 



coincident avec les developpements des resultats exacts obtenus par d'autres methodes 
pour le front de diffusion |1141 1221 , pour la largeur thermique et pour le parametre 
de localisation |3^. Ces comparaisons avec des resultats exacts obtenus par d'autres 
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methodes montre que la procedure de renormalisation generalisee est exacte ordre par 
ordre en /x. Plus generalement, pour calculer les observables a I'ordre //", il suffit de 
considerer que le front de diffusion s'etale sur au plus (1 + n) pieges et de moyenner 
sur les echantillons avec la mesure appropriee. Et comme nous avons une description 
du front de diffusion echantillon par echantillon, cette methode permet de calculer 
perturbativement toutes les observables que Ton souhaite (Publication [Pll]). 

Cette approche permet par ailleurs de bien comprendre comment, dans la phase de 
diffusion anormale < fi < 1, ilyaala fois une largeur thermique qui croit en i^'^ (|4.4|1 
et une probabilite finie l2(/u) (|4.5|) de trouver deux particules au meme point a temps 
infini. Notre description de la repartition du paquet thermique dans un echantillon 
donne est tout a fait en accord avec les simulations numeriques de A. Compte et J. P. 
Bouchaud ^^, dont la figure reproduite ici H4.2() montre que meme pour fi = 0.4, il 
n'y a typiquement que quatre pieges qui jouent un role important dans un echantillon 
donne a un instant donne. 

4.4 Equivalence quantitative entre le modele de Sinai bi- 
aise et le modele dirige de pieges 



Fig. 4.3 - Calcul du temps de sortie d'une vallee renormalisee de barriere T : V ex- 
pression exacte du temps de premier passage en b pour une particule qui part de fait 
intervenir une integrate double qui est dominee par le facteur d'Arrhenius e^^ , avec un 
prefacteur qui fait intervenir deux fonctions de partition : Zy represente la fonction 
de partition de la vallee et Zb represente la fonction de partition du potentiel (—1^) 
autour du sommet de la barriere P. 

Une analyse par la methode du col du temps de sortie d'une vallee renormalisee 
donnee de barriere P montre qu'il est distribue selon une loi exponentielle, comme dans 
le modele de piege, avec un temps de piege 

Q ~ ^ZsZye^^ (4.6) 
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qui depend surtout de la barriere F a travers de facteur d'Arrhenius usuel e^^, mais qui 
depend aussi des details de la vallee a travers le prefacteur qui contient deux fonctions 
de partition de vallees Browniennes independantes Zy et Zb f Figure 113)) . 

La distribution du temps de piege Ty sur I'ensemble des vallees du paysage renor- 
malise a I'echelle F pent etre calculee a partir de la distribution des barrieres et du 
prefacteur : on retrouve alors la meme expression que dans le paysage renornialise du 
modele dirige de piege 

q,(r) = 9{t < r)f^ l^iy (4.7) 

a condition de choisir I'echelle de renormalisation F du paysage de Sinai en fonction du 
temps selon 

r {t) = T In \ta^/3^{r^ {I + fj.))^] (4.8) 

L'echelle de longueur associee 

m = ^[ta'pr (4-9) 

correspond alors exactement a la constante du front de diffusion qui a ete calculee 
recemment par les mathematiciens probabilistes j93j . 

Finalement, notre approche etablit que le modele de Sinai biaise et le modele 
dirige de pieges sont asymptotiquement equivalents du point de vue de leurs descrip- 
tions renormalisees a grande echelle, a une echelle de longueur pres que nous avons 
determinee. 

4.5 Resultats pour le modele de Sinai biaise 

Tons les resultats du modele dirige de pieges peuvent done etre traduits pour le 
modele de Sinai biaise, en remplagant les pieges ponctuels par des vallees renormalisees. 
Le front de diffusion dans un echantillon fixe a une structure hierarchique analogue, 
representee par la Figure 14.41 qui est I'equivalent de la figure 14.11 Dans les resultats 
quantitatifs, il suffit de remplacer la variable d'echelle X = ^ par la variable d'echelle 

X = |W adaptee au modele de Sinai avec biais (|4.9jl . En particulier, la largeur ther- 
mique admet le developpement suivant 



< Ax2(t) > (^2^3) 



2/o3\2m 



t2M cj2/34 



(2 In 2) , 7r2 , ,, ^1^1 

^^+ -— + 21n2 1n2-2-27s K + - 



(4.10) 



Et si on traduit le resultat exact obtenu pour le modele dirige de pieges [B], on obtient 
meme le resultat suivant pour toute la phase < // < 1 de diffusion anormale 

< Ax2(t) > _ (a2/?3)2^r4(^) ^sinvr/.^^^^^ ^^^^^^ 



t2M (j2;g4 r(2^) V ^Ai 
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Marche aleatoire dans un potentiel Brownien biaise 




Fig. 4.4 - Structure hierarchique des vallees importantes pour une particule partant 
de I'origine. Les barrieres opposees au biais, qui sont soulignees par des lignes droites, 
correspondent aux profondeurs du modele de pieges de la Figure \4 1\ Le fond M de la 
vallee renormalisee qui contient I'origine a I'echelle T est avec une probabilite d'ordre 
0(1). Le fond Li de la vallee renormalisee suivante et le fond Si de la plus grande sous- 
vallee avant M sont occupies avec des poids d'ordre 0{fj,). Le fond L2 de la deuxieme 
vallee renormalisee apres M , la plus grande sous-vallee I2 entre M et Li, et la seconde 
plus grande sous-vallee S2 avant M sont occupees avec des probabilites d'ordre 0(/x^). 



en ternies de I'integrale 



Hi^) 



^2^+2 _^ 2 COS 7^|J.z^^+^ + 1 



(4.12) 



4.6 Conclusion 

La phase de diffusion anormale x ~ t^ avec < /i < 1 dans le modele de Sinai 
biaise est caracterisee par une localisation sur plusieurs vallees renornialisees, dont les 
positions et les poids peuvent etre caracterises dans chaque echantillon. La procedure 
de renormalisation generalisee permet de calculer toutes les observables souhaitees par 
un developpement perturbatif exact en fi. 

Du point de vue des methodes de renormalisation de type Ma-Dasgupta, cette etude 
montre que la procedure usuelle qui est exacte lorsque le desordre evolue vers un point 
fixe de desordre infini, est une approximation qui garde un grand interet lorsque le point 
fixe est caracterise par un desordre fini assez grand : la procedure usuelle constitue alors 
I'ordre dominant d'un developpement perturbatif systematique. 
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Chapitre 5 

Chaines de spin classiques 
desordonnees 

5.1 Presentation des modeles 

L'equilibre thermodynamique des chaines de spins classiques desordonnees peut etre 
formule comme un produit de matrices de transfert 2x2 aleatoires. En particulier, 1' 
energie libre par spin correspond a I'exposant de Lyapunov et peut done etre etudiee 
par la methode de Dyson- Schmidt : le livre de J.M. Luck |123j expose les nombreux 
resultats obtenus dans ce cadre. Dans ce chapitre, nous allons etudier l'equilibre et la 
dynamique de ces modeles, en termes des domaines Imry-Ma. 

5.1.1 La chaine d'Ising en champ aleatoire 

La chaine d'Ising en champ aleatoire a pour Hamiltonien 

n=N-l i=N 

n = -J Y, SiSi+i-Y,hiSi (5.1) 

Les champs {hi} sont des variables aleatoires independantes, de moyenne nulle hi = 
(cf la publication [P7] pour le cas d'une moyenne non-nulle ) et de variance 

g^hl (5.2) 

D'apres I'argument de Imry-Ma |lU7j . I'etat fondamental a temperature nulle est desordonne 
il y a une alternance de domaines de spins (+) et de domaines de spins (— ), la taille 
typique d'un domaine etant la longueur de Imry-Ma 

LiM ~ (5.3) 

9 

L'argument de Imry-Ma est le suivant jlU7j : la creation d'un domaine de taille L coute 
une energie 4 J (creation de deux parois), qui est independante de L, mais elle permet de 
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gagner une energie typique |2^^^^ ^i|typ ~ 2^^^ . En consequence, pour L > Ljm, il 
devient energetiquement favorable de creer un doniaine pour profiter d'une fluctuation 
favorable des champs aleatoires. 

5.1.2 La chaine verre de spin en champ magnetique exterieur 

L'Hamiltonien du verre de spin en champ magnetique s'ecrit 

i=N-l N 

W = - ^ JiCTiCTi+i - ^ hai (5.4) 

ou les couplages { Jj} sont des variables aleatoires independantes. Dans le cas particulier 
Ji = lb J avec probabilite (1/2, 1/2), il existe une equivalence avec la chaine en champ 
aleatoire (J5.1J1 avec /ij = it/i : 

- On pose Ji = Jei avec ej = ±1 

- On effectue une transformation de jauge des spins ai = ei-.e-i-iSi 

- On pose hi = /iei..ej_i 

L'interpretation physique de cette correspondance entre les deux modeles est la 
suivante : 

• En champ nul h = 0, les deux etats fondamentaux du verre de spin zbCj- , avec 
(Tj. = ibei..ej_i, correspondent aux etats ferromagnetiques de la chaine d'Ising pure 
5(0) = +1 et 5f ) = -1. 

• En presence de /i > 0, I'etat fondamental est une succession des deux fondamen- 
taux de champ nul de taille typique Ljm = ^^■ 

• Les parois de domaines de la chaine en champ aleatoire representent les liens 
frustres JiUiUiJ^i = JSiSi+i < du verres de spin. 

Dans la suite, nous ne decrirons done les resultats que dans le langage de chaine 
en champ aleatoire, car il est immediat de les traduire pour la chaine verre de spin en 
champ exterieur 

5.2 Principe de la renormalisation 

5.2.1 Hamiltonien pour les parois de domaines 

L'Hamiltonien (|5.1() pour les spins correspond a I'Hamiltonien suivant pour les 
parois de domaines ^o(+|— ) et Ba{—\+) dans I'ordre AiBiA2B2... : 

Na Nb 

n = 2J{NA + NB) + Y.V{ac:)-Y,V{ho) (5.5) 

a=l a=l 

en terme du potentiel de Sinai 

X 

V{x) = -2^hi (5.6) 

En consequence : 
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- Chaque parol (A ou B) coute une energie 2 J 

- Les parois Aa(+|— ) voient le potentiel V{x) 

- Les parois Ba{—\+) voient le potentiel oppose {—V{x)) 

5.2.2 Dynamique de croissances de domaines 

On s'interesse a la dynamique de Glauber a partir d'une configuration initiale 
aleatoire : le taux de transition est W{Sj — > —Sj] 



— e/5AB_|_g-/3AE avec 
AE{ creation de 2 parois} = 4 J it 2hj 
AE{ diffusion d'une parol} = ziz2hj 
AE{ annihilation de 2 parois} = — 4J it 2hj 



(5.7) 



Dans le regime {hi} ^ T <^ J qui va nous interesser dans toute la suite de ce 
chapitre, la croissance de domaines depuis / ~ 1 jusqu'a Ljm ^ 1 est decrite par le 
modele suivant de Reaction-diffusion dans le potentiel de Sinai 

- Les parois A diffusent vers les minima 

- Les parois B diffusent vers les maxima 

- Annihilation lors des rencontre A + i? — > 

On pent de plus montrer (Publication [P7]) qu'a grand temps, tons les maxima et 
minima sont occupes par des parois, comme sur la Figure ()5.1() . ce qui simplifie I'analyse 
du processus de reaction-diffusion par la renormalisation du paysage. 




A 



Fig. 5.1 - La ligne en zig-zag represente la renormalisation du potentiel de Sinai i5. 6]) 
vu par les parois \5.5\) . Les parois de domaines de type A{+\—) occupent les minima, 
alors que les parois de domaines de type B{—\-\-) occupent les maxima. Les descentes 
contiennent done des spins (+), alors que les montees contiennent des spins (— ). 



5.2.3 Equilibre 

Le processus de reaction-diffusion decrivant la croissance de domaines s'arrete a 
I'echelle de renormalisation 



Peg = Tin teq = 4 J 



(5.^ 
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lorsqu'on atteint 1' equilibre avec la creation de paires de domaines. 

5.3 Resultats 

5.3.1 Densite de parois 



i5_ 
(rint)2 



La densite de parois decroit selon n(t) = ^rr^„ ^x-?. jusqu'a I'equilibre Ueg = 1/Ljm- 



5.3.2 Distribution des longueurs des domaines 

Ce sont des variables independantes (ce qui n'est pas vrai pour le cas pur 63 
distribuees selon la loi 



CO / 1 \ 2 

P*(A)=7r Y, (n+-j(-l)"e-"'^("+5) 



(5.9) 



avec la variable d'echelle A = /J,^^^^^ au cours de la croissance de domaines, et A = ^ — 
a I'equilibre. 

5.3.3 Fonction de correlation spatiale 



^ 48 + 64(2n + I)27r2c/M .,„ , -,^2^2oJ-l 

{So{t)sAt)) = Y. — d^;TT)i^^' ^^-^^^ 

n=~oo ^ ' 

avec r = (Tint) au cours de la croissance de domaines, et Teg = 4 J a I'equilibre. 

5.3.4 Fonction de correlation temporelle 

La fonction d'autocorrelation d'un spin decroit selon 

4 /Inf'A 1 /Int'A^ 



ce qui correspond a un exposant d'autocorrelation a grand temps A = 1/2. 

5.3.5 Exposants de persistance 

• La probabilite qu'un spin ne se retourne pas dans [0, t] a une decroissance 

((Tw) avec0 = L 

• La probabilite que la valeur moyenne thermique < Si(t) > ne change pas de signe 
dans [0, t] decroit selon : 

3-V5 



La probabilite qu'un domaine initial n'ait pas disparu a I'instant t decroit selon : 

{Tlnty\ 



' avec ip = y . (De maniere generale, on a I'inegalite ip < 
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5.3.6 Violation du theoreme de fluctuation-dissipation 

5.3.6.1 Rapport X pour mesurer le caractere "hors equilibre" 

Le parametre X de violation du theoreme de Fluctuation-Dissipation est defini par 

m 

T R{t,t^)=Xit,t^)dt^Cit,t^) (5.12) 

oil C{t,tw) represente la fonction de correlation thermique tronquee 

Cit,t^) = Y,<So{t)S^{t^) >-< So{t) X S,{t^) > (5.13) 

X 

et oil R{t, tu]) represente la fonction de reponse lineaire lorsqu'on applique un champ 
uniforme H a partir de t^ 

<So{t)> = H I duR{t,u) (5.14) 

Pour la chaine en champ aleatoire, on trouve trois regimes 
(i) un quasi-equilibre des parois dans les vallees 

X{t,t^) = l pour < ^"^,^ ~ ^"^^ < 1 (5.15) 

(ii) un rapport X non trivial lorsque la dynamique effective des vallees repart 

X{t,tw) = — pour — fixe (5.16) 

(iii) un regime final de vieillissement 

twlntw 7 In t i-ntu, 

En particulier, X croit vers +oo, car les correlations thermiques tronquees sont tres 
faibles par rapport a la reponse a un champ. 

5.3.6.2 Comparaison avec les mo deles de champ moyen 

Ici, le rapport X n'est pas une fonction de la correlation C{t,tw) : c'est a cause des 
deux echelles {t,tw) et (lnt,lnt^„)). 

D'autre part, en champ moyen le rapport X appartient a I'intervalle [0, 1] et s'in- 
terprete comme I'inverse d'une temperature effective X = 1/Tf.ff |15] . 

Ici, on trouve T^ff — > 0, ce que Ton pent interpreter comme un point fixe de 
temperature nulle. 
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5.3.6.3 Comparaison avec la croissance dans les ferromagnetiques purs 

II est interessant de comparer avec ce qui est connu sur le regime de croissance dans 
les ferromagnetiques purs : 

• Au point critique T = Tc, il existe un rapport X{-^) non-trivial qui correspond a 
un rapport d'amplitude ^2| Par exemple, pour la chaine d'Ising pure unidimensionnelle 
a temperature nulle Tc = 0, le rapport d'amplitude |H3] vaut X{t, t^) = ^-^ : il decroit 
de X{t = tw,tw) = 1 a X{t — > oo, tw) = g- 

• Pour T < Tc, le rapport X est nul X = ^1 E], ce qui est interprete de 
la maniere suivante : les parois repondent par un facteur 0(1) mais elles n'occupent 
qu'une fraction 1/L(t^) du volume. 

Par comparaison, on voit que pour la chaine en champ aleatoire, seule une petite 
fraction l/P^ des parois repondent, mais avec une reponse tres grande, qui correspond 
au renversement de tout un domaine d'ordre r?„. 



w 



5.4 Conclusion 

La formulation des chaines de spins classiques desordonnees en termes de parois 
de domaines qui voient un potentiel de Sinai, permet d'etudier a la fois les proprietes 
d'equilibre thermodynamique et la dynamique de croissances de doinaines a partir d'une 
condition initiale aleatoire. La possibilite d'obtenir des resultats exacts vient de la tres 
forte localisation des parois de domaines par le desordre. 

Publications associees : 

• Principe de la renormalisation et premiers resultats : Publication [P3] 

• Etude detaillee : Publication [P7] 

• Etude de processus de reaction-diffusion plus generaux dans un potentiel de Sinai : 
Publication [P5]. 



Chapitre 6 

Localisation d'un polymere 
aleatoire a une interface 

6.1 Presentation du modele 

Ce chapitre est consacre a un modele introduit par Garel, Huse, Leibler et Orland 
[SJ : un heteropolymere constitue de monomeres de charges qi aleatoires, est en presence 
d'une interface situee en z = qui separe un miheu z > favorable aux charges positives 
g > et un milieu z < favorable aux charges negatives q < 0. Plus precisement, la 
version continue de ce modele est definie par la fonction de partition suivante sur les 
trajectoires Browniennes {z{s)} j81j 

ZL{P;{q{s)})=Jvz{s)exp(-^l^ ds (^^^ + P J^ dsq{s)sgn{z{s))] (6.1) 

A haute temperature, des arguments de type Imry-Ma ont ete proposes |S2 pour le cas 
symetrique gj = et pour le cas biaise gj > : ce sont des arguments energie/entropie, 
a la difference des arguments usuels d'Imry-Ma qui sont de type energie/ energie. 

6.1.1 Argument de type Imry-Ma pour le cas symetrique 

L'argument de type Imry-Ma pour le cas symetrique [SJI est le suivant. On suppose 
que la chaine est localisee autour de I'interface, avec des boucles de longueur typique / 
dans chaque solvant : 

i+l 

• L' energie typique gagnee par boucle est d'ordre y^ qi ~ vol 

i 

• La perte d'entropie pour une boucle est d'ordre (In /) 
L 'optimisation de 1' energie libre par monomere 
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par rapport a la longueur / conduit a 

I T2 



(6.3) 



Get argument predit done une localisation du polymere symetrique a toute temperature, 
avec les comportements d'echelle suivants pour la longueur typique d'une boucle et 1' 
energie libre 

l(T) ~ ^ ^ (6.4) 

a 

f{T) ~ --^ (6.5) 

6.1.2 Argument de type Imry-Ma pour le cas dissymetrique 

L'argument de type Imry-Ma pour le cas dissymetrique [H^ est en fait beaucoup 
plus subtil que dans le cas symetrique, parce que la description correcte des boucles 
dans le solvant (— ) necessite la consideration des "evenements rares" oii la somme 
de I- variables aleatoires qi, de moyenne positive gj = go > 0, se trouve etre assez 
negative pour rendre favorable une excursion dans le solvant (— ). Plus precisement, 
l'argument propose [SJ est le suivant : on s' attend a ce que le polymere soit dans le 
solvant prefere (+), sauf lorsqu'une boucle de longueur /~ dans le solvant (— ) devient 
favorable energetiquement, avec une charge Q- = — Ylil=j Qi > ^ suffisante. Comme 
la probabilite d'avoir J2l=j Qi = ~Q 

Prob(g") = e i^F (6.6) 

est faible, la distance typique l^ entre deux tels evenements se comporte comme I'inverse 
de cette probabilite 

1+ n^ e 4ai- . (6.7) 

Get argument probabiliste de type 'evenement rare', qui est inhabituel dans les argu- 
ments de type Imry-Ma qui concernent les evenements typiques, donne que 1' energie 
gagnee Q~ dans une boucle du solvant (— ) se comporte en 

Q- r^ Vial- In 1+ -qor , (6.8) 

alors que la perte d'entropie associee est d'ordre (In l^). La difference d' energie libre par 
monomere entre cet etat localise avec des boucles (/+, Z_) par rapport a I'etat delocalise 
dans le solvant prefere (+) est d'ordre 

f{T,l+,l-)-fdeioc{T) ~ p (-Q- + rin/+) ~ p (gor-V4cTMn/+ + rin/+) (6.9) 
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L 'optimisation par rapport a la longueur l~ donne 

r~-^lnZ+ (6.10) 

% 

et done finalement 

Q-~— ln/+ (6.11) 

Ainsi, cet argument predit qu'a la fois 1' energie gagnee Q~ et le cout en entropie ont 
la meme dependance en (Inl^) : la difference d' energie libre se factorise selon 

f{T, /+) - fdeiociT) ~ (T - re)i^ (6.12) 

Cet argument conduit done finalement a une transition a la temperature critique 

T, ~ - (6.13) 

entre la phase localisee T < Tc et la phase delocalisee T > T^.- Contrairement au cas 
symetrique, les comportements de 1' energie libre et de la longueur /+ par rapport a 
la temperature ne sont pas determines par I'argument. L'objet de la renornialisation 
de type Ma-Dasgupta va etre justement de trouver ces comportements en construisant 
explicitement les boucles d'Imry-Ma du polymere autour de I'interface en fonction de 
la realisation du desordre. 

6.2 Principe de la renornialisation 

6.2.1 Construction de la structure optimale en boucles 

A temperature nulle T = 0, chaque monomere veut etre dans son solvant prefere 
sgn(2;j) = sgn((7j) : le polymere est decompose en boucles a contenant la monomeres 
consecutifs de meme signe, et portant une certaine charge absolue Qa- Lorsque la 
temperature T croit, on considere les configurations de la chaine qui peuvent etre 
obtenues a partir de la structure de I'etat fondamental en transferant dans le solvant 
oppose les boucles de plus petite charge absolue Qmin = T de maniere iterative. Quand 
on transfere une boucle (Q2 = Tj h) qui est entouree de deux boucles {Qi, h) et (Qs, ^3), 
on obtient une nouvelle boucle de charge absolue Q et de longueur I donnees par les 
regies (cf Figure 16. Ij) 

Q = Qi + Qi-Q2 (6.14) 

l = li + h + h 

Ces regies correspondent aux regies de renornialisation des extrema Browniens 
decrites dans le Chapitre [21 : ici la marche aleatoire correspondante est la somme des 
charges ^q qj en fonction du monomere final i, et la procedure de renormalisation 
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(QA) 



z=0 



(Q,,l,) 



^ 



(Qyl,) 







z=0 




Fig. 6.1 - Illustration de la procedure de renormalisation dans I'espace reel : le transfert 
d'une boucle {Q2 = ^jh) entouree par les boucles voisines (QiJi) et {Q3J3) donne 
naissance a une nouvelle boucle {Q,l) avec les regies de renormalisation if). 15]) . 



construit la meilleure structure en boucles avec la contrainte qu'on n'autorise que les 
boucles de charge absolue plus grande que T. 

Pour etablir la correspondance entre I'echelle de renormalisation T et la temperature 
T, il reste a etudier dans quelles conditions le transfert d'une boucle ((52; '2) dans le 
solvant oppose est effectivement favorable pour 1' energie libre : le coiit en energie est 

AEf^'P = 2Q2 (6.15) 

alors que le gain en entropie est 

ASf^'P = IniMih +I2 + k)) - ln[M{li)Mil2)M{l3)] (6.16) 

ou Ai{l) = S /l^'"^ represente le nombre de marches aleatoires de / pas qui vont de z = 
a z = en presence d'une frontiere absorbante en z 
energie libre 



, ce qui conduit au bilan d' 



^pihp ^ ^E^UP _ TAS^^'P = 2Q2 - Tin 



Mjh + h + h) 
M{h)M{l2)M{h) 



(6.17) 



Pour obtenir la structure optimale a temperature T, il faut done poursuivre la renormal- 
isation tant qu'elle permet d'abaisser 1' energie libre [AF^^'^p < 0), et il faut s'arreter 
a la premiere iteration qui conduirait a une augmentation (AF-' ^ > 0). On definit 
Teq{T) comme I'echelle de renormalisation T oil il faut arreter la renormalisation. 



Resultats pour le cas symetrique 53 

6.2.2 Validite de la methode 

La methode de renormalisation est justifiee si I'echelle de renormalisation Teq{T) 
est grande, ce qui correspond a la region des hautes temperatures. 



6.3 Resultats pour le cas symetrique 

6.3.1 Energie libre 

L' energie libre trouvee 

fiT) ~ -^ (6.18) 

est en accord avec I'argument de Imry-Ma ()6.5p . 



6.3.2 Distribution des longueurs de boucles 

Les longueurs de boucles sont des variables aleatoires independantes, et la variable 
d'echelle 

(6.19) 



9r2(lnT)2 
est distribuee avec la loi 



/ 1 \ 2 2 

P(A) = Yj (n + -)vr(-l)"e-^'^("+^) ^~ Tve-^^ (6.20) 



A— >oo 



^ E (-i)'"(™ + 5''"*"'"*'\:„^^"* (^^^" 



rra=— oo 



6.3.3 Densite de polymere autour de I'interface 

La variable d'echelle pour la distance z a I'interface 




D 3T{\nT) 
est distribuee avec la loi 



(6.22) 



R{Z) = 4 / dXP{X)Vx / due'''' ~ -^2^6^'^^ (6.23) 

Jo J^ ^^°° ^ 

V A 
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6.3.4 Temperature de delocalisation d'une chaine finie 

Alors qu'une chaine infinie reste localisee a toute temperature, une chaine finie 
de taille L a une dissymetrie d'ordre yL pour sa charge totale et subit done une 
delocalisation a une temperature finie. Pour caracteriser ces effets de taille finie, on 
pent calculer la distribution de la temperature de delocalisation T^eioc sur I'ensemble 
des chaines cycliques de taille finie grande L : la variable d'echelle 

3 

g = —F^Tdcioc In Tdeioc (6.24) 

VcrL 

a pour distribution 

D(g) =^> (-l)"+^n2e "iP" ~ — e ^ (6.25) 

= ^ y r2(2m + l)V-l]e-(2"+^)'^' ^ ^ffV^' (6.26) 

V rfr^ . — r^v^ y 



6.4 Resultats pour le cas dissymetrique 

6.4.1 Temperature critique 

Dans le cas dissymetrique (go > 0), dans la limite a ^ qo, la temperature critique 
obtenue par renormalisation 

Tc = ^ (6.27) 

est dans le domaine des hautes temperatures, et done dans le domaine de validite de la 
methode. 

6.4.2 Singularite essentielle de 1' energie libre 

La renormalisation conduit a la singularite essentielle suivante pour 1' energie libre 



In^ 

go 



1- ^ 



(6.28) 



f{T)-f{T,) ~ _2gofln— )exp 
T->T- V 90 / 

ce qui correspond done a une transition d'ordre infinie. 

6.4.3 Proprietes des boucles dans les deux solvants 

La longueur typique l^igi,{T) des boucles dans le solvant prefere diverge avec une 
singularite essentielle a la transition, alors que la longueur typique l^i^^ijiT) des boucles 
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dans le solvant defavorable diverge algebriquement 



a 
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90 



T^T- %' f 1 - J 



(6.29) 
(6.30) 



La variable d'echelle A+ = /^//^^^^(T) pour les boucles dans le solvant prefere est 
distribuee avec la distribution exponentielle e~ + . 



6.4.4 Densite de polymere autour de I'interface 

La variable d'echelle pour la distance z a I'interface dans le solvant (+) 

7. 

z = 



f^^-p 






a pour distribution 

R+{Z) = A / 
Jo 
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(6.31) 



(6.32) 



6.4.5 Temperature de delocalisation d'une chaine finie 

Pour caracteriser les effets de taille finie, on pent calculer la distribution de la 
temperature de delocalisation Tdeioc sur I'ensemble des chaines cycliques de taille L : la 
variable aleatoire 



d / 4cr\ Tc-Tj;,io^ 



QqL V Qo 



est distribuee selon la loi 



D+( 



1 _i 

—pre r 



(6.33) 



(6.34) 



En particulier, la valeur typique pour la temperature de delocalisation d'une chaine 
de taille L presente une correction logarithmique d'ordre (1/lnL) par rapport a la 
temperature critique Tc de la limite thermodynamique 
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6.5 Conclusion 

La renormalisation de type Ma-Dasgupta permet d'etudier la localisation d'un 
heteropolyniere aleatoire autour d'une interface a haute temperature, en construisant 
explicitement la succession des boucles dans chaque echantillon. Ces boucles sont des 
domaines d'Imry-Ma generalises de type energie/entropie. Dans le cas dissymetrique, 
la competition entre les energies aleatoires gagnees dans les boucles et les entropies per- 
dues lors des retours a I'interface conduit a une transition de delocalisation caracterisee 
par des singularites essentielles. 
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Chapitre 7 

Potentiel Brownien en presence 
d'un confinement quadratique 

7.1 Presentation du modele 

Le potentiel aleatoire unidimensionnel 

W-) = f.' + VW (T.l) 

contenant un terme quadratique deterministe et un terme aleatoire Brownien V{x) 



{V{x)-V{y)f = 2\x-y\ (7.2) 

a ete introduit par Villain et al. comme un 'modele jouet' pour les interfaces en presence 
de champ aleatoire |163j . Par ailleurs, dans le cadre de I'etude des varietes de dimension 
interne D plongees dans un milieu aleatoire de dimension (A^ + D), il est considere 
comme le cas extreme le plus simple D = et N = 1 |134j . 

Ce modele simple permet en effet de mieux comprendre certaines proprietes de 
systemes desordonnes plus compliques. 

7.1.1 Argument de Imry-Ma 

A cause du confinement quadratique, le minimum absolu du potentiel est fini et 
un argument de Imry-Ma, entre 1' energie elastique d'ordre ^x^ et 1' energie aleatoire 
d'ordre ^/x, permet de determiner son ordre de grandeur 



-"mm 



/x-2/3 (7.3) 



alors que les methodes usuelles de perturbation a tons les ordres jl64j ou d'iteration 
jl65j n'arrivent pas du tout a reproduire ce scaling ()7.3() . Dans le cadre de la methode 
variationnelle des repliques, le resultat H7.3|) necessite une brisure de symetrie des 
repliques pHllTHj. 
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7.1.2 Symetrie statistique de 'tilt' 

Ce modele presente aussi une symetrie statistique de 'tilt' ( 'statistical tilt sym- 
metry' en anglais) comme les autres modeles en champs aleatoires, qui implique des 
identites remarquables |154j pour les moyennes thermiques des cumulants thermiques 
de la position, que Ton pent resumer par jl54j 



A 



2 



In < e-^^- > = T — (7.4) 

Cette identite sur la fonction generatrice montre que le second cumulant est simplement 



<x2>-<x>2 = — (7.5) 

et que toutes les moyennes sur le desordre des cumulants superieurs s'annulent ! Dans 
le regime de basse temperature, le resultat (|7.5|) implique que les fluctuations ther- 
miques sont reliees a la presence d'etats metastables dans des echantillons rares |154j . 
La renormalisation permet en particulier d'etudier quantitativement ce phenomene. 



7.2 Principe de la renormalisation 

Pour un potentiel aleatoire U{x) 'Markovien', satisfaisant une equation de Langevin 
locale de type 

^^ = F[C/(x),:E]+r?(x) (7.6) 

ax 

ou r/(x) est un bruit blanc Gaussien, la mesure du paysage renormalise se factorise en 
blocs qui satisfont des equations de renormalisation fermees (Publication [P9]). Pour le 
cas de paysages stationnaires, oil la force F est independante de x 

m.l = FPl = -^ (7.7) 

(ce qui generalise le paysage Brownien pur F = et le paysage Brownien biaise F[U] = 
F > 0), et pour le modele (|7.1j) qui correspond a une force independante de U et lineaire 
en X 

F[U{x),x]=F[x]=nx (7.8) 

on obtient des solutions explicites pour la renormalisation (Publication [P9]). 

En particulier, pour le modele (|7.1|) qui nous interesse dans ce chapitre, la mesure 
du paysage renormalise s'exprime en terme des fonctions d'Airy. 
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7.3 Resultats sur la statistique des minima 

7.3.1 Position du minimum 

A temperature nulle, la particule se trouve au minimum Xmin du potentiel. L'etat 
final r = cxD de la renormalisation permet de retrouver que la distribution de Xmin sur 
I'ensemble des echantillons est, en accord avec PSIIHO] 

P\-oo,+ool{x) = g{x)g{-x) (7.9) 

et la fonction auxiliaire 

en utilisant les notations a = (/x/2)"'^''^ et 6 = 1/a^, 

7.3.2 Echantillons avec deux minima presque degeneres 

La probabilite qu'un echantillon presente deux minima presque degeneres avec une 
difference d' energies AE = e ^ situes en xi et X2 s'ecrit 

'D{e,xi,X2) = eg{-xi)d{x2 - xi)g{x2) +0{e^) (7.11) 

en termes de la fonction g ()7.1U|) et de la fonction 

La probabilite d'avoir deux minima separes par une distance y > est done 

En particulier, le calcul du second moment donne 

r+oo 1 

/ dyy^D{y) = - (7.14) 

Jo /^ 

7.3.3 Contribution au second cumulant a basse temperature 

La contribution des echantillons avec deux minima degeneres au second cumulant 
thermique de la position (|7.5() pent etre estime a I'ordre T en temperature comme suit : 
les deux minima ont pour poids de Boltzmann respectifs p = ^ L^^ and {1 — p) = 

^^^-0e ■ La variable {x— < x >) est done (1 — p)(xi — X2) avec probabilite p et p(x2 — xi) 
avec probabilite 1 — p, et done apres moyenne on obtient 



< {x- < X >)'^ > = p{l - p){xi - X2y (7.15) 

/+00 r+00 r+00 e~''l'^ 

dx\ \ dx2 \ (ieP'(e = 0,xi,X2)— — — ?yr^(x2 -xi)^ (7.16) 
-00 J xi J —00 (i + e ) 



(1+6-^/^)2 

T I dyD{y)y' (7.17) 

/o 



' —00 J X\ 
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En utilisant H7.14() . on obtient alors exactement le resultat exact 1)7. 5|) . Ceci montre 
que les fluctuations thermiques a basse temperature viennent entierement des etats 
metastables qui existent dans quelques echantillons rares. En particulier, la suscepti- 
bilite 

X = ^ (< x2 > - < X >2) (7.18) 

a une moyenne finie a temperature nulle 

X=- (7.19) 

mais seuls les echantillons avec deux minima degeneres contribuent a cette valeur 
moyenne, car les echantillons typiques avec un seul minimum ont une susceptibilite 
qui s'annule a temperature nulle. 

7.3.4 Contribution aux moments thermiques a basse temperature 

De meme, les moments pairs de la position relative (x— < x >) se comportent de 
la maniere suivante a basse temperature 



T 



+00 



<{x-<x >)2« > = - / y'^'Diy) + 0{T') (7.20) 

n Jo 

en terme de la fonction D(y) definie en H7.13() . La comparaison avec I'identite ()7.4|) 
montre qu'il y a beaucoup de termes qui se compensent dans les moyennes sur le 
desordre des cumulants. 



7.4 Resultats sur la statistique de la plus grande barriere 

Le temps igg necessaire pour atteindre I'equilibre est directement relie a la plus 
grande barriere Tmax = Tlntgg qui existe dans I'echantillon. Plus precisement, la prob- 
abilite Viteq < t) que le systeme ait deja atteint I'equilibre a I'instant t, correspond a 
la probabilite qu'il ne reste plus qu'une vallee renormalisee a I'echelle F = Tint : c'est 
une fonction de la variable d'echelle 7 

nte, < t) = $ (^7 ^ (I) '^' Tint) (7.21) 

qui a pour expression explicite en termes des fonctions d'Airy 

.(,) ^/;;^^e-/o -/..(/.) (7.22) 

J,(fX) = Ai{f + j + iX) 

^^^■'' ^ 7rAi{f + iX)[Ai{f + iX)Bi{f + -f + iX)-Bi{f + i\)Ai{f + -f + iX)] 
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Ce resultat correspond directement par derivation a la distribution de probabilite 
de la variable reduite 7 = (^) Tmax de la plus grande barriere Tmax de Pechantillon 

Pmaxh) = ^^h) (7.23) 

Les comportements asymptotiques de cette distribution sont les suivants 

Pmaxh) - Ixl^l'^'e-h'^' (7.24) 

7^00 4 V 2 

PmaAl) =^,^e '^ (7.25) 

Oil C(^) designe la fonction zeta de Riemann. 

7.5 Conclusion 

Pour le potentiel Brownien avec confinement quadratique, la solution explicite du 
paysage renormalise en termes de fonctions d'Airy permet d'etudier la statistique des 
minima et des barrieres. En particulier, ceci permet de montrer que les fluctuations 
thermiques a basse temperature sont entierement gouvernes par les echantillons rares 
qui possedent deux minima degeneres. 

Publication associee [P9] 
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Chapitre 8 

Le modele de pieges 
unidimensionnel 

8.1 Presentation des modeles de pieges 

Les modeles de pieges proposent un mecanisme tres simple de vieillissement [2S| '■ 
une particule effectue une dynamique stochastique dans un paysage contenant des pieges 
d' energies aleatoires de distribution exponentielle 

piE)=e{E)^e-^ (8.1) 

-'a 

Ce choix de distribution exponentielle est justifie par les etats de basse energie du REM 
[S5] . la theorie des repliques jl33j . et plus generalement par la queue exponentielle de 
la distribution de Gumbel qui represente une classe d'universalite importante dans la 
statistique des extremes. Cette distribution exponentielle des energies (|8.1() correspond 
pour le temps de piegeage d'Arrhenius r = e^^ a la loi algebrique 



qir)=9ir>l)-^^ (8.2) 



avec I'exposant 



l^ = 7fr (8.3) 

-'-9 

A basse temperature T < T^, le temps moyen de piegeage J dTTq{T) diverge, ce qui 
conduit immediatement a des effets de vieillissement. 

Les proprietes de vieillissement ont ete beaucoup etudiees, soit dans la version de 
champ moyen [^Sl El HH I11U| - soit dans la version unidimensionnelle [751 ^\, qui 
apparait naturellement dans divers contextes physiques j^ H21 (HHl ; et qui presente deux 
echelles de temps caracteristiques pour le vieillissement. L'objet de ce chapitre est de 
comprendre ce phenomene grace a une renormalisation appropriee. 
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Le modele de pieges unidimensionnel 



8.2 Principe de la renormalisation 

Nous avons deja decrit la renormalisation adaptee au modele dirigede pieges, dans le 
Chapitre0J Ici, dans le modele non-dirige , chaque site pent etre visite un grand nombre 
de fois, ce qui conduit a un changement essentiel : un piege du paysage renormalise va 
etre caracterise par deux temps importants, a savoir 

(i) son temps de piegeage Tj, qui represente le temps typique de sortie vers ses voisins 
immediats 

(ii) son temps d'evasion Tj, qui represente le temps necessaire pour arriver dans un 
piege plus profond. 



Fig. 8.1 - Definition du temps d'evasion dans le paysage renormalise : le piege de 
temps de piegeage tq est entoure par deux pieges, le piege r+ a la distance /+ et le piege 
T- a la distance /_. Le temps d'evasion Tq est le temps caracteristique necessaire pour 
atteindre r+ on r_ quand on part de tq. 



Plus precisement, le paysage renormalise a I'echelle R est defini comme suit : tons 
les pieges Ti < R sont remplaces par un paysage 'plat', alors tons les pieges Ti > R 
sont inchanges. Dans le paysage renormalise (cf Figure IH?T|l . quand la particule quitte 
le piege tq, elle s'echappe vers la droite ou vers la gauche avec probabilite 1/2. Si elle 
s'echappe vers la gauche, elle va arriver a atteindre le piege suivant r_ avec probabilite 
l//_. Si elle s'echappe vers la droite, elle va arriver a atteindre le piege suivant t+ avec 
probabilite 1/1+. Sinon, elle va etre reabsorbee par le piege tq. Asymptotiquement, le 
nombre de retours dans le piege tq avant I'evasion va done etre grand, d'ordre R^. On 
pent montrer que le temps passe a I'interieur de tq pendant ces multiples visites est 
dominant par rapport au temps passe dans les multiples excursions infructueuses et 
dans I'excursion reussie pour s'echapper. Finalement, on obtient que le temps d'evasion 
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tgsc du piege tq a une distribution exponentielle 



1 -Ih 



P{tesc)^^ T^e ^0 (8.4) 



avec le temps caracteristique 



i+ + L 

ce qui explique pourquoi il y a deux echelles de temps importantes. En particulier, 
pour decrire la dynamique au temps t, on souhaite garder les pieges de temps d'evasion 
T > t dont la particule n'a pas eu le temps de s'echapper, et eliminer au contraire 
les pieges de temps d'evasion T < t, ce qui conduit au choix suivant pour I'echelle de 
renormalisation R{t) du paysage en fonction de temps 

R{t) ~ tW (8.6) 

De plus, dans la limite fi ^ 0, la dynamique effective suivante devient exacte : a 
I'instant t, la particule partie de I'origine se trouve soit sur le premier piege renormalise 
M+ a distance /+ sur la droite, soit sur le premier piege renormalise M_ a distance l- 
sur la gauche. Le poids de M+ est simplement la probabilite 1-/(1^ + 1-) d'avoir atteint 
M+ avant M_ 

Peffix, t) ~ j^^Hx + 1-) + J^^i^ - h) (8-7) 

8.3 Resultats 

8.3.1 Distance entre pieges dans le paysage renormalise 

Les distances entre pieges successifs dans le paysage renormalise sont des variables 
aleatoires independantes et la variable reduite A = l/S,{t) est distribuee selon la loi 
exponentielle 

P(A) = e-^ (8.8) 

L'echelle de longueur caracteristique du paysage renormalise a I'instant t est 

m = 6(/^) t^ (8.9) 

avec le prefacteur 

eo(/x) = l + 0(/x) (8.10) 
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8.3.2 Forme d'echelle du front de diffusion 

Dans la variable d'echelle X = x/^{t), la moyenne sur les echantillons du front de 
diffusion (|8?7|) est 

r+co 

g^{X) = e-l^l / due-"^-—— + 0(/i) (8.11) 

Jo |A| + tt 

8.3.3 Parametres de localisation 

Les parametres de localisation, qui representent les moyennes sur les echantillons 
des probabilites que k particules independantes se trouvent au meme site, sont donnes 
par 



+00 2 



Ykif,) ^ hm J2 ^'(«' i|0, 0) = --— - + 0(/x) (8.12) 

t-+oo ^-^_ {k + 1) 



n=0 



ce qui est en accord avec les simulations numeriques de Bertin et Bouchaud ^7j qui 
ont obtenu Y2 — > 2/3 et Y3 -^ 1/2 dans la limite ^ ^ 0. 

8.3.4 Fonction generatrice des cumulants thermiques 

La largeur thermique est 



< n^ > — < n >^ 



C2(^)=lim -27- = l + 0{fi) (8.13) 

et plus generalement, les autres cumulants thermiques peuvent etre derives de la fonc- 
tion generatrice 

jj f+00 / sA sA \ 

Z^(s) = In < e""«w > = / dAe^^Af — coth— -1] +0(//) (8.14) 

8.3.5 Fonction de correlation a deux particules 

La fonction de correlation a deux particules prend la forme 



+00 +CO ^ / 7 \ 



n=0 m=0 

Le poids de la fonction 6 a I'origine correspond comme il se doit au parametre de 
localisation Y2 = 2/3+0(/i) (|8.12|) . alors que la seconde partie est une fonction d'echelle 
de la variable A = t4y qui a pour expression 

C^(A) = e-^^ + 0(/i) (8.15) 
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8.3.6 Les deux fonctions de correlations temporelles 

La probabilite n(t + t^, tw) de ne pas sauter pendant I'intervalle de temps [t«,, t«, + t] 
a une forme d'echelle de 'sous-vieillissement' 

t \ ^ f ,^ , ,,_L t 

r 



n(t + 1», t^) = n^ ( <7 = ;^^7^ ) = n^ I 5 = [To(/i)] 1+^ ^ ) (s.ie) 



avec 



n(0)(^) = /" dz^z^-ig-^s (8.17) 

En particulier, on obtient le comportement asymptotique 

n(t + t^,t^) ~ (^) [/^ + 0(/ti^)] (8.18) 

— ^^ >+oo \ , 1+p / 

,1+71 ^^w ' 

'"W 

La probabilite C{t + t-u,^tw) d'etre au temps (t + t«,) dans le piege occupe a I'instant 
tw a une forme d'echelle de vieillissement 

C{t + tw,t^) = C^(h = j^r^^) =C,(^h = To(/i)^) (8.19) 

avec 

/2h 



Cf\h) = C,{h) = j^j^ dzz'+^'^Kfiz) (8.20) 

En particulier, on obtient le comportement asymptotique 

C{t + tw,tJ ~ (^) [^ + 0(/w')] (8.21) 

T^ >oo \tw / 

Dans cette approche de renormalisation, la presence de deux echelles de temps 
differentes dans les proprietes de vieillissement vient directement de la presence d'un 
temps d'evasion different du temps de piegeage. Comme en dimension d = 2 une marche 
aleatoire est encore recurrente, on s'attend la aussi a avoir deux echelles de temps, 
alors que pour d > 2, il n'y a plus de recurrence, et la seule echelle de temps pour le 
vieillissement est I'echelle des temps de piegeages. II est done interessant de discuter le 
cas d = 2. 

8.4 Analyse de scaling pour le modele de pieges en li = 2 

8.4.1 Paysage renormalise 

On definit encore le paysage renormalise a I'echelle R en ne gardant que les pieges 
de temps de piegeage Tj > R. Etant donne un piege, la distance / au plus proche voisin 
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a pour distribution 

27r/ 

Pr{1)^^ -j^e R^ (8.22) 

En utilisant le fait que la probabilite de s'echapper jusqu'a une distance / sans etre 
revenu a I'origine se comporte en c/(lnZ), on obtient que le nombre n de retours avant 
une evasion reussie se comporte en n ~ In/ '^ lni?2'. Finalenient, le temps d'evasion 
d'un piege tq du paysage renormalise est de la forme 

ro = arolni?t (8.23) 

ou a est une variable aleatoire d'ordre 1 qui caracterise la geometric des pieges voisins 
de Tq. L'echelle de renormalisation R doit done etre choisie en fonction du temps selon 

tr^R{t)\nR2{t) (8.24) 

ce qui donne par inversion a I'ordre dominant 

R{t) ~ -^ (8.25) 

Inta 

8.4.2 Deux echelles de temps 

La presence de deux echelles de temps distinctes R{t) et t va conduire a deux 
comportements de vieillissement : on s'attend a ce que la correlation n(t + tw,tw) soit 
une fonction de t/R{tw) avec 1)8. 25() et a ce que la correlationC(t + tw,tw) soit une 
fonction de t/tw, ce qui est en accord avec les resultats rigoureux fT2| [T3]. 

8.4.3 Deux echelles de longueur 

La grande nouveaute qualitative par rapport au cas d = 1 est qu'il existe deux 
echelles de longueur en d = 2. 

8.4.3.1 Distance entre deux pieges renormalises 

La distance caracteristique entre deux pieges voisins du paysage renormalise est 



m-[Riw 



t 



lnt^' 



A» 



(8.26) 



8.4.3.2 Distance parcourue a I'instant t 

Par ailleurs, I'analyse de scaling a partir des sommes de Levy J2J en d = 2 est la 
suivante : apres A^ pas, le nombre de sites distincts visites croit selon S = N/{lnN) et 
chaque site a ete visite typiquement (In A^) fois, ce qui conduit a la correspondance 

t r^J2^,{lnN) ~ (IniV) (^) " (8.27) 
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Comme la distance typique apres N pas est r ~ v A'^, on obtient finalement le com- 
portement en temps de la distance parcourue 



1- 



r(t) ~t2(lnt'^)^" (8.28) 

8.4.3.3 Nombre de pieges importants a I'instant t 

En d = 2, la presence des deux echelles de longueur differentes ^(t) (|8.26|) et r{t) 
H8.28() et ^(t) ()8.26() montre que le nombre n{t) de pieges renormalises contenus dans le 
disque de rayon r{t) accessible, ne reste pas fini mais croit logarithmiquement en temps 

"W-^-l"*'' (8-29) 

ce qui est en accord avec le resultat rigoureux concernant I'absence de localisation en 
d = 2 a temps infini t ^^ oo JSj. Cependant, dans le regime t — > oo et /x ^^ avec 
t^ fixe H8.29() . on s'attend a ce que le front de diffusion dans un echantillon soit une 
somme de pics 6 dont les poids dependent de la repartition geometrique des pieges mais 
pas de leurs energies. II serait bien sur interessant d'arriver a definir une dynamique 
effective plus precise pour decrire ce regime, a partir d'une construction explicite des 
probabilites d'absorption des differents pieges pour une configuration aleatoire fixee. 

8.5 Conclusion 

La renormalisation appropriee pour decrire la phase de vieillissement // < 1 du 
modele de pieges unidimensionnel est basee sur deux echelles de temps, le temps de 
piegeage et le temps d'evasion, qui gouvernent deux fonctions de correlations tem- 
porelles differentes. A I'ordre dominant en ^u ^ 0, dans un echantillon donne, le front 
de diffusion est localise seulement sur deux sites, qui correspondent aux deux pieges 
renormalises les plus proches de I'origine a I'echelle de renormalisation R{t). Les poids 
des pieges sont simplement donnes par les probabilites respectives d'atteindre un piege 
avant I'autre. Ceci montre que la dynamique reste toujours hors equilibre : les poids 
des deux pieges ne sont pas donnes par des facteurs de Boltzmann, ils ne dependent 
meme pas des energies des pieges, mais seulement des distances a I'origine. 

Publication associee [P12] 
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Chapitre 9 

Le modele de pieges en champ : 
reponses lineaire et non-lineaire 

9.1 Presentation du modele 

En presence d'un champ exterieur /, I'equation maitresse du modele de pieges s'ecrit 

m 

'-^=Py\n + l)^+Py\n-l)^-Py\nf^^l±^ (9.1) 
at 2t„+i 2r„_i 27:„ 

dans laquelle les taux de transition satisfont la condition de bilan detaille 

-/3t/(n)^{/) -/3C/„+ip^{/) (9_2) 

en terme de 1' energie totale C/„ contenant 1' energie aleatoire {—En) du piege et 1' 
energie potentielle lineaire {—fn) associee au champ exterieur / 

Un = -En - fn (9.3) 

La reponse de ce modele a ete etudiee par E. Bertin et J. P. Bouchaud ^Hl avec 
des arguments de scaling et des simulations numeriques. Leur resultat principal est 
le suivant : le Theoreme Fluctuation-Dissipation du regime de reponse lineaire est 
valable meme dans le secteur de vieillissement, mais la reponse devient non-lineaire 
asymptotiquement. Le but du present chapitre est de comprendre ces proprietes, en 
utilisant d'une part une symetrie particuliere de I'equation maitresse, et d'autre part 
en generalisant la renormalisation du chapitre precedent pour prendre en compte la 
presence d'un champ exterieur. 

9.2 Une symetrie dynamique qui empeche la violation du 
theoreme Fluctuation-Dissipation 

Dans cette section, les equations seront ecrites pour le cas d= \ pour simplifier les 
notations, mais il est clair que les resultats peuvent etre generalises immediatement a 
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toute dimension finie d. 

9.2.1 Loi de conservation dans chaque echantillon pour le niodele 
sans champ 

Le modele de pieges a une propriete tres speciale US] : en I'absence de champ 
exterieur / = la moyenne thermique de la position 

+ 00 

<n>f=o{t)^ Y, nPy='\n) (9.4) 

71= — OO 

est une quantite conservee dans tout echantillon ! Cela vient du fait qu'en sortant d'un 
piege la particule saute avec des probabilites egales (1/2, 1/2) vers la droite ou vers la 
gauche. Cette propriete de conservation est d'autant plus remarquable que le desordre 
brise la symetrie n — > —n dans chaque echantillon. 

9.2.2 Un theoreme de fluctuation non-lineaire pour t^ = 

Dans un echantillon donne, si on compare les fronts de diffusion en presence de (+/) 
et de (— /), on obtient la relation tres simple 



i+f)i 

e^^" (9.5) 



pr\n)_^,f^ 



Pt^\n 
ce qui permet d'exprimer la difference des moyennes thermiques de la position 

+00 

< n >+/ (t)- < n >„/ (t) 
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Pt^\n) 



(9.6) 



9.2.2.1 Relation de Fluctuation-Dissipation dans le regime lineaire 

La premiere consequence ()9.6|) est que la relation d'Einstein 

<n>f (t) ^^^ ^<n^ >/=o (t) (9.7) 

est valable pour t fixe dans la limite / ^ 0. 

9.2.2.2 Fin du regime de reponse lineaire 

La deuxieme consequence de (|9.6p est que pour / fixe, le regime de reponse lineaire 
n'est valable que s'il est legitime de lineariser le facteur {l — e~^^"') ce qui donne comme 
critere 

/5/e(/=o)(i)«l (9.8) 



Une symetrie dynamique qui empeche la violation du theoreme Fluctuation-Dissipation 73 

ou C{f=o)i't) est le deplacement typique sans champ a Pinstant t. Le regime de reponse 
lineaire est done limite par un temps caracteristique associe a la force / 

l + M 
1 



t « W) ^ [jj) (9.9) 

9.2.2.3 Propriete du front de diffusion 

Une troisieme consequence de ()9.6() concerne la dissymetrie de la moyenne sur les 
echantillons du front de diffusion (Publication [P13]) 



Py\-n) = e-'^f''Py\n) (9.10) 

9.2.3 Un theoreme de fluctuation non-lineaire pour t^ quelconque 

Les resultats precedents peuvent etre generalises au cas oii le systeme evolue d'abord 
sans champ exterieur / = pendant I'intervalle [0,t^], avant d'etre soumis au champ 
/ pour t > tw On introduit le front de diffusion a deux temps P^^''^^'{n,t]nw,tw\^,Q) 
qui represente la probabilite jointe d'etre en n^y a i^, et en n a i (on suppose t > t^)- 

La generalisation de H9.5() 

p(+/-*-)(n,t;n^,t^|0,0) ^ ^/3/(„_„„) ,g ^^. 

P(-/;*-)(n,t;n^,t^|0,0) 

permet d'obtenir les trois consequences suivantes. 

9.2.3.1 Regime lineaire pour (t, i«,) fixes 

Pour des temps (t,t«,) fixes, la relation d'Einstein est valide sous la forme 

0f 
<{n- n^) >f,t^ (t) ~ ^ < (n - n^f >/=o (0 (9.12) 

9.2.3.2 Fin du regime de reponse lineaire 

Pour un temps t^ fixe et un champ exterieur / fixe, la validite du regime lineaire 
est limitee aux temps t tels que 

/3/e/=o(t,t«,)«l (9.13) 

ou (,f=o{t,tw) represente le deplacement caracteristique entre t^ et t en I'absence de 
champ, et done finalement 

t-t^^t^{f) (9.14) 

ee qui generalise H9.9() . 
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9.2.3.3 Front de diffusion a deux temps 

L'asymetrie induite par / sur du front diffusion moyen est simplenient 



P(/)(-n,t;-n^,t^|0,0) =e-'3^("-"™)p(/)(n,t;n^,t»|0,0) (9.15) 

9.2.4 Discussion 

Si on considere une equation maitresse generale, la condition pour deriver un theoreme 
de fluctuation non-lineaire analogue a celui qui existe pour le modele de pieges est la 
suivante (Publication [P13]) : pour toute configuration C, le taux de transition global 
pour en sortir en presence de (+/) doit etre egal au taux de transition global pour en 
sortir en presence de (— /) 

Wg^l (C) = w^^l (C) pour toute configuration C (9.16) 

A I'ordre lineaire en /, cette condition implique que I'observable b lineairement couple 
au champ / (qui generalise la position dans le modele de pieges) satisfait une loi de 
conservation : sa moyenne thermique doit etre conservee en I'absence de champ pour 
une condition initiale arbitraire. 

Existe-t-il d'autres modeles interessants qui possedent cette symetrie ? La fagon la 
plus simple de satisfaire la condition (|9.1()j) serait de pouvoir associer a toute transition 
C ^ C une autre transition C -^ C" qui a le meme taux de transition en champ nul 
W^^'{C — > C") = W^^'{C -^ C) et telle que la variation de b soit exactement I'opposee 
h{C") — b{C) = —{b{C') — b{C)). Cependant, c'est une propriete tres forte car elle doit 
etre vraie pour toute configuration initiale C. 

9.3 Etude de la reponse par renormalisation 

Dans le Chapitre precedent sur le modele de pieges sans biais, nous avons decrit une 
procedure de renormalisation a partir de la notion du temps d'evasion et des probabilites 
respectives de s'echapper vers le piege renormalise de droite ou vers le piege renormalise 
de gauche. En generalisant cette approche a la presence d'un champ exterieur /, on 
pent obtenir beaucoup de resultats exacts sur la dynamique. 

9.3.1 Resultats explicites pour t^ = ^ 

La renormalisation permet de definir une echelle de longueur ^(i, /) qui represente 
la distance moyenne entre deux pieges dans le paysage renormalise associe a I'instant 
t. Cette echelle interpole entre les deux comportements 

^(i,/) :^ tTl^[l + 0(/.)] (9.17) 

t«t^(/) 

e(t,/) :^ (/3/t)qi + 0(/i)] (9.18) 
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9.3.1.1 Front de diffusion 

La moyenne sur les echantillons du front de diffusion prend la forme d'eclielle 

1 



PtinJ) 



;9, X 



n 



-,F = f5mt,f) 



et la fonction d'echelle g^ a pour expression, dans la limite /i ^ 0, 



(9.19) 



go{X, F) = e-l^l [^(X > 0) + 9{X < 0)e-^l^l] 
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-\X\ 





+ 00 



+°° _;, 1 - e- 

1 - e-^(l^l+^) 



9{X > 0) + 9{X < 0)e-^l^l] Y, Jl 

n=0 ^ 



F 



(9.20) 

-nF|X| 



{I + Fn){l + F + Fn) 



9.3.1.2 Position moyenne 



La moyenne sur les echantillons de la moyenne tiiermique de la position est de la 
forme 



+ 00 



<x{tj)>^ Y, xPtixj) = at,f)MF = pmtj)) 

X = — OD 

et la fonction d'echelle X^ a pour expression, dans la limite /x ^ 0, 

f + OO 

MF) -- 



(9.21) 



+".^^ A^COth^-1 
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^--.-^*"^--. 



(9.22) 



9.3.1.3 Largeur thermique 

La moyenne sur les echantillons de la largeur thermique est de la forme 

<Ax2(t,/)> = fit, /)A^(F = (3mt, /)) 
et la fonction d'echelle Aq est en fait directement reliee a la fonction Xq 



MF) 
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MF) 



(9.23) 



(9.24) 



9.3.2 Resultats explicites pour tt„ > 

9.3.2.1 Dynamique effective en fonction de {f,t,tu,) 

L'etat atteint au temps t^ en I'absence de champ doit etre considere comme une 
condition initiale pour la dynamique en presence du champ / > en fonction du 
nouveau temps (t — tw). Comme a I'instant t^, la particule est typiquement dans un 
piege r > R{tw, / = 0), la dynamique effective ne recommence que lorsque la procedure 
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de decimation redevient active pour R{t — tw,f) > R{tw,f = 0). II est done utile 
d'introduire le parametre 

qui mesure les echelles de longueur des deux paysages renormalises a t^; et a t. Dans 
le secteur a < 1, la reponse est domine par des evenements rares, alors que dans le 
secteur a{t,tyj,f) > 1, la reponse est domine par la dynamique effective. Le domaine 
en {t — tiu) correspondant au secteur a > 1 depend des valeurs relatives de t^ et du 
temps t^{f). 

(i) Pour t^ <C t^{f), le secteur a > 1 correspond au domaine temporel (t — tw) > i^, 
et le parametre a se comporte selon 

(t — t \ 1 + M 
^^ for t„<t-tw«t^(f) (9.26) 

et 

/il(t-t )\'" 
a{t,t^J) = 2 v^^; for t-tw>t^(f) (9.27) 

1 
(ii) Pour t^ » ^m(/)' ^^ secteur a > 1 correspond au domaine {t — t^) > -wjtlu^^ , et 

le parametre a se comporte comme H9.27|l partout. 

En etudiant la statistique du paysage renormalise a deux echelles de renormalisation 

successives, on pent alors obtenir des resultats explicites valables dans tout le secteur 

a > 1 dans la limite /i ^ 

9.3.2.2 Loi du deplacement entre t^ et t 

En termes des variables d'echelle Y = ^}u_^ J\ , F = PfS,{t — t«,5 /) et a = ^{t — 
tw, f)/C{tw, f = 0), la moyenne sur les echantillons de la distribution de Y est de la 
forme 

P(Y;F,a) = -5{Y)+ \e{Y >0) + e-^\^\e{Y < 0)] Gns{\Y\,F,a) (9.28) 

La partie singuliere en 5 represente les particules qui n'arrivent pas a s'echapper vers 
un autre piege renormalise entre ty^ et t. La partie reguliere est definie en termes de la 
fonction 

Gns (Y, F,a) = ^- J^ ^^ ^_^_^(y+„) [(2 - e-("-^)^)e-- - e"""] (9.29) 

Plus generalement, on pent calculer ainsi toutes les fonctions d'echelle souhaitees 
dans la limite /i — > 0, en particulier pour la position moyenne et la largeur thermique 
(Publication [P14]). 
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9.4 Conclusion 

C'est la symetrie particuliere de I'equation maitresse du modele de pieges qui 
empeche la violation du Theoreme de Fluctuation-Dissipation nieme dans le regime 
de vieillissement. Par ailleurs, la renormalisation basee sur les temps d'evasion, qui de- 
vient exacte dans la limite /x ^ 0, permet d'etudier en detail la reponse non-lineaire 
dans les differents regimes qui existent en fonction des deux temps (tyj , t) et du champ 
exterieur /. 

Publications associees 

Sur le theoreme de Fluctuation-Dissipation non-lineaire : [P13] 
Sur les calculs explicites par renormalisation : [P14] 
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Chapitre 10 

Chaine de spin quantique S = 1 
antiferromagnetique aleatoire 

10.1 La chaine antiferromagnetique S* = 1 sans desordre 

10.1.1 Differences entre spin demi-entier et spin entier 

Alors que la chaine antiferromagnetique S = 1/2 presente des correlations en loi 
de puissance et des excitations sans gap, la chaine antiferromagnetique 5 = 1 est car- 
acterisee par des correlations a decroissance exponentielle et par un gap pour les exci- 
tations. Une fagon simple de comprendre ces differences qui existent plus generalement 
entre spins demi-entiers et spins entiers est la fonction d'onde "Valence-Bond-Solid" 
(VBS), qui permet de mettre en evidence un ordre a longue portee pour un parametre 
d'ordre topologique qui est non-local en termes de spins. 

10.1.2 La fonction d'onde VBS 

Si on represente chaque spin S = 1 comme la symetrisation de deux spins elementaires 
S = 1/2, la fonction d'onde VBS [J est I'etat obtenu en format un singulet sur chaque 
lien entre deux spins 1/2 constitutifs (Figure riU.l|l . Cette fonction d'onde est I'etat 



n-n-n-Ti-rt 



Fig. 10.1 - Fonction d'onde VBS : chaque spin 5 = 1 est represente par la symetrisation 
de 2 spins 1/2 (les boules noires reliees par un trait vertical gras) ; les lignes horizontales 
qui relient deux spins 1/2 representent les singulets entre spins 1/2 constitutifs. 
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fondamental exact de I'Hamiltonien pP 

HaKLT = Yl ^2(^, + S^+l) = Y, 
i i 

ou P2 est le projecteur sur le sous-espace s = 2 



-Si.Si+i + -{Si.Si+i) +- 
lb 6 



flO.l) 



10.1.3 Le parametre d'ordre topologique 

Le parametre d'ordre |145j qui permet de caracteriser I'ordre topologique de la 
fonction d'onde VBS est non local dans les spins 



Uj = -(V'olSf exp 



m 



Y.SI 



i<k<j 



^fko) 



flO.2) 



-N-^ 



Pour I'etat pur VBS sur une chaine cyclique de N spins, on obtient tij = 4/9 + 0(3" 
Pour le fondamental de la chaine antiferromagnetique pure, ce parametre d'ordre tend 
aussi vers une valeur finie, ce qui indique une parente avec la fonction d'onde VBS. 



10.2 Construction des regies de renormalisation 

10.2.1 Renormalisation d'un lien antiferromagnetique 

L' Hamiltonien d'un lien antiferromagnetique entre deux spins S = 1 



hn 



J1S1.S2 — — 



^\2 

5*1 + 5*2 



S' 



c2 
'-'2 



Jl 

2 



S1 + S2 



presente trois niveaux d' energie indexes par la valeur s = 0, 1, 2 du spin total 



J^ 



[sis + l)-i] 



(10.3) 



(10.4) 



il y a done un singulet eo = — 2Ji, un triplet ei = — Ji et un quintuplet 62 = Ji- 

Si on generalise naiVement la regie de Ma-Dasgupta du cas S = 1/2 en projetant 
sur le singulet, le nouveau couplage effectif est 



J^ 



4J0J2 



3 Jl 



(10.5) 



Le coefficient 3 etant plus grand que 1, cette regie n'est pas automatiquement consis- 
tante, et la procedure ne pent etre justifiee que si on part d'un desordre assez large. 

Nous avons done propose de generaliser la procedure de Ma-Dasgupta avec la regie 
suivante : au lieu de projeter sur le niveau le plus bas de ho, il faut en fait eliminer 
le niveau le plus haut. Ainsi, pour I'hamiltonien ho, on elimine le quintuplet, mais 
on garde le singulet et le triplet, en remplagant les deux spins S = 1 par deux spins 
S = 1/2. Cette decimation partielle elargit I'espace initial, mais il est possible de definir 
une renormalisation fermee avec 4 types de liens. 



Construction des regies de renormalisation 



81 



10.2.2 Procedure de renormalisation avec 4 types de liens 

La procedure de renormalisation que nous avons proposee est definie pour 1' ensem- 
ble elargi des chaines constituees de spins de taille S = 1/2 ou 5 = 1, dans lesquelles 
les couplages {Jj} sont soit ferromagnetiques (F) soit antiferromagnetiques (AF), avec 
la contrainte suivante : pour tout segment {i,j}, 1' aimantation classique doit satisfaire 
I^TiijI < 1. Cette condition pour deux voisins j = i + 1 montre qu'il y a 4 types de liens 
possibles : 

1) Lien F entre deux spins S=l/2, 

2) Lien AF entre deux spins S=l/2, 

3) Lien AF entre un spin S=l et un spin S=l/2, 

4) Lien AF entre deux spins S=l. 

Les quatre regies de renormalisation correspondantes sont les suivantes : 



(1) 
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So Si = i S2 = ^ S3 
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So Si = 1 S2 = ^ 
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10.2.3 Interpretation de la renormalisation en termes d'amas VBS 

Si on represente chaque spin 5=1 initial comme la symetrisation de deux spins 
S = 1/2, les regies 2, 3, et 4 pour les liens AF s' interpretent comme la formation d' un 
singulet entre deux spins 5 = 1/2 constitutifs. La regie 1 pour un lien F entre deux spins 
S = 1/2 correspond a leur symetrisation. A la fin de la procedure de renormalisation, 
lorsqu il n'y a plus de spin libre, la chaine est decomposee en un ensemble d'amas 
disjoints qui ont une structure VBS ( Figure [Tn.2|) . 

Pour la chaine desordonnee, on a done un parametre d'ordre (|lfl.2|) t^ = 4/9 si les 
deux sites appartiennent au meme amas et tij = sinon. Pour une chaine de taille N, 
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Chaine de spin quantique 5 = 1 antiferromagnetique aleatoire 




Fig. 10.2 - L'etat fondamental obtenu par renormalisation, en termes d'amas VBS : 
chaque spin S = 1 est represente par la symetrisation de 2 spins 1/2 (les boules noires 
reliees par un trait vertical gras) ; les lignes qui relient deux spins 1/2 representent des 
singulets. Les spins 1 initiaux sont ainsi regroupes en amas : ici, par exemple, il y a 
deux amas de taille 2, un amas de taille 4, et un amas plus grand qui depasse les limites 
de la figure. 

la moyenne spatiale ^^ • tij/N"^ du parametre tij est proportionnelle a la probabilite T 
que deux spins appartiennent au meme amas 



r-E|l = !4E'«+o(i/iV) 



(10.6) 



*J 



Oil Uc est le nombre de spins dans un amas c, qui n'est pas directement liee a son 
extension spatiale. Ce parametre d'ordre T ne pent etre non-nul dans la limite thermo- 
dynamique que s'il existe un amas VBS qui contient une fraction finie des spins de la 
chaine. 



10.3 Etude numerique de la renormalisation 

Nous avons etudie numeriquement la procedure de renormalisation avec 4 types 
de liens en partant de chaines cycliques de N spins (par exemple A^ = 2^^ ~ 4.10^ ), 
avec des couplages initiaux Jj uniformement distribues dans I'intervalle [1,1 + d]. Le 
parametre d represente done la largeur du desordre initial. Pour chaque taille, nous 
avons effectue des moyennes sur un certain nombre d'echantillons (typiquement 100) . 
Nous avons etudie le Hot des quantites suivantes en fonction de 1' echelle F : 
(i) le nombre A^(F) de spins effectifs S = 1/2 et S = 1 encore presents a 1' echelle F 
(ii) la proportion {N (^g^i-^(j--^ / N iV)] de spins S = 1 parmi les spins effectifs presents 
(iii) les proportions piiV) = {A'^j(F)/A^(F)} des liens de type i = 1,2,3,4 
(iv) les distributions de probabilite Pi{J.,Q) des couplages J pour les quatre types 
de liens i = 1,2, 3,4. 

Les resultats de la renormalisation presentent un changement qualitatif pour une 
certaine valeur critique dc — 5.75(5) du desordre initial. Ainsi, la Figure ITO .31 represente 

le flot de la proportion '■^/p\ — de spins 5 = 1 pour differentes valeurs du desordre 
initial : il y a deux valeurs attractives, a savoir a faible desordre et 1 a fort desordre. 



Proportion of spin-1 among the effective spins 

for initial disorder d=1 , 2, 4, 5.5, 6, 8,16,1 00 




50.0 



Fig. 10.3 - Proportion de spins S = 1 parmi les spins ejjectifs a I'echelle T pour 
differentes valeurs du desordre initial d = 1,2,3,4,5.5,6,8,16,100 ; cette proportion 
tend vers dans la phase de faihle desordre, et vers 1 dans la phase de fort desordre. 
Au point critique dc — 5.75(5) la proportion reste stationnaire autour de la valeur 
intermediaire 0.315(5). 



10.3.1 Resultats dans la phase de fort desordre 

Dans la phase de fort desordre d > do le nombre A^(r) de spins effectifs decroit 
comme dans la phase appelee "Random Singlet" de la chaine S = 1/2 : 



'^('^'r^^F- 



;io.7) 



et les proportions PiiT) des quatre types de liens convergent vers un regime asympto- 
tique caracterise par les proportions simples (Fig 110.41) 

Pi(r)~0 p2(r)~e(r) p^{T)^2e{T) p4(r) ~ 1 - 3e(r) (10.8) 

oil e(r) tend lentement vers en F — > oo. II y a done presque partout des liens de 
type 4, avec quelques defauts de type (lien de type 3,lien de type 2, lien de type 3) qui 
viennent de la decimation partielle momentanee des liens de type 4. A fort desordre, la 
renormalisation converge done vers la phase "Random Singlet" de Ma-Dasgupta. 



10.3.2 Resultats dans la phase de faible desordre 

Dans la phase de faible desordre, le nombre A^(F) de de spins effectifs decroit expo- 
nentiellement 



NiV) oc e 
r^oo 



-o(d)r 



(10.9) 



avec un coefficient a{d) qui est une fonction decroissante de d qui s' annule a la transition 
d ^> d~ . Les proportions piiT) des quatre types de liens convergent vers le regime 
asymptotique (Figure ri().5|l 



Pi(F)~o.25 p2(r)~o.75 P3(r)~o P4(r)~o 



(10.10) 
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Chaine de spin quantique S = 1 antiferroniagnetique aleatoire 

Proportions p,{r), p^{r), p^iT)/!, p^{r) 

for large initial disorder (d=100) 
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Fig. 10.4 - Les proportions Pi{T) des 4 types de liens i = 1,2,3,4 of en fonction de 
Vechelle T pour un desordre initial fort d = 100. 
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Proportions p,{T), p^CH, P3(r), p^cn 

for weak initial disorder (d=0.1) 
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Fig. 10.5 - Les proportions pi{T) des 4 types de liens i = 1,2,3,4 of en fonction de 
Vechelle T pour un desordre initial faihle d = 0.1 
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10.3.3 Resultats numeriques au point critique 

10.3.4 Etude des types de liens 

Au point critique, le nombre de spins efFectifs decroit algebriquement 

r— >oo i •^ 



(10.11) 



et les proportions Pi{T) des quatre types de lien convergent vers le regime asymptotique 
(Figure IMH) 



pi(r)~0.17, 



P2(r) ~0.35, 



P3(r)~0.33. 



P4(r)~0.15. (10.12) 
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Fig. 10.6 - Les proportions Pi{T) des 4 types de liens i = 1,2,3,4 of en fonction de 
Vechelle T pour un desordre initial critique dc = 5.75 . 



10.3.5 Etude numerique de la transition de percolation des amas 
VBS 

Du point de vue des amas VBS, la transition de phase quantique correspond a 
une transition de percolation : dans la phase de fort desordre, il n'y a que des amas 
finis, alors que dans la phase de faible desordre, il apparait un amas macroscopique qui 
contient une fraction finie des spins. 



10.3.5.1 Parametre d'ordre 

Si on note (3 I'exposant qui gouverne 1' annulation de la fraction ni/N des spins 
dans I'amas macroscopique, le parametre d'ordre se comporte en T ~ (dc — d)"^^ pour 
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d < dc- L'etude de finite size scaling de la Figure ()1U.7|) conduit a I'estimation 

2/3 = 1.0(1). (10.13) 
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Fig. 10.7 - Parametre d'ordre topologique T en fonction du desordre initial d pour des 
chaines de tailles N = 2^^ — 2^°. Le point critique est situe en dc = 5.76(2) 



10.3.5.2 Susceptibilite 

La Figure H10.8|) represente la taille moyenne des amas finis, qui joue le role d'une 
susceptibilite 



v = V^ 



Ol 



(10.14) 



(c=l est I'amas le plus grand). Si on note 7 I'exposant critique qui gouverne sa diver- 
gence X ~ \dc — d\~'^ , on obtient I'estimation 



7 = 1.2(1). 



(10.15) 



10.3.5.3 Distribution des amas au point critique 

Au point critique, la distribution mc{s) de la taille s des amas a une decroissance 
algebrique fFigure riU.9|) 



VflrXs) 



(10.16) 



avec un exposant 



r = 2.2(1). 



(10.17) 



Calcul des exposants critiques exacts 
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Fig. 10.8 - Susceptibilite x e?^ fonction du desordre initial d pour des chaines de tailles 

AT = 2^4 -220. 




10 10 10 

size s 



Fig. 10.9 - Distribution rnds) de la taille s des amas an point critique : la mesure de 
la pente conduit d Vestimation r = 2.2(1) 



10.4 Calcul des exposants critiques exacts 

II existe un modele effectif plus simple de chaine de spin jHE] qui se trouve etre 
dans la meme classe d'universalite (Publications [PI]). En introduisant une variable 
auxiliaire /i pour etudier la taille des amas VBS dans ce modele plus simple, nous 
avons calcule les exposants exacts de la transition 



2;3 = lt^ = 1.17278... 



7 



\/l3-l 
2(2^/5-3) 
^/13-1 



1.13000..., 



r = 1 + ^ = 2.34164... 

V5 



(10.18) 

(10.19) 
(10.20) 



qui sont en tres bon accord avec nos estimations numeriques ()10.13p . H10.15|) . 1)10.1711 . 
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10.5 Conclusion 

La generalisation de la niethode de Ma-Dasgupta a la chaine de spin S = 1 antifer- 
romagnetique aleatoire conduit a un modele a 4 types de liens que nous avons etudie 
numeriquement. La transition de phase observee en fonction du desordre initial s'in- 
terprete comme une transition de percolation pour les amas VBS de I'etat fondamental. 
En utilisant un modele efFectif dans la meme classe d'universalite, nous avons calcule 
les exposants critiques exacts de cette transition de phase quantique. 

Publications associees 

Sur la transition de percolation des amas VBS : [PI] 

Sur les degres de liberte efFectifs de basse energie dans les differentes phases : [P2] 



Chapitre 11 

Chaine d'Ising avec couplages et 
champs transverses aleatoires 

11.1 Presentation du modele 

11.1.1 Transition de phase quantique 

La chaine d'Ising quantique 

H = -Y^ JiSfSf^^ - Y, h,Sf (11.1) 

i i 

avec couplages Jj > et champs transverses hi > aleatoires (les signes peuvent etre 
fixes par une transformation de jauge), est le modele desordonne le plus simple qui 
presente une transition de phase quantique a temperature nulle : le point critique situe 
en 



ln/i = lnJ (11.2) 



separe une phase ferromagnetique (In/i < In J) d'une phase desordonnee (In/i > In J). 
Ce modele pent etre etudie en detail par la renormalisation de type Ma-Dasgupta [TSj 
comme nous I'avons deja explique dans I'lntroduction (Page|^. Dans ce Chapitre, nous 
allons nous interesser aux proprietes de taille finie au voisinage de la transition, dans 
deux ensembles possibles pour le desordre. 

11.1.2 Ensembles microcanonique et canonique pour un systeme 
desordonne fini 

Dans I'etude des systemes desordonnes, on considere generalement que les variables 
aleatoires de desordre dans un echantillon sont des variables independantes : cet ensem- 
ble d'echantillons sera appele ici 'canonique' comme dans les references jnHIMl- (ce qui 
correspond a I'ensemble appele 'grand canonique' dans d'autres references |1471 El 11681 
[T3]). Cependant, il est apparu qu'il pouvait etre interessant jl47j de considerer aussi 

89 
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un ensemble d'echantillons appele ' micro canonique' dans les references |991 I66j ( qui 
correspond a I'ensemble appele 'canonique' dans les references |147| ini ri68| ll5j). Dans 
cet ensemble microcanonique, il existe une contrainte globale sur les variables aleatoires 
de desordre dans un echantillon de A^ sites. Le premier exemple important concerne un 
systeme pur avec une fraction p £ [0, 1] d'impuretes |147j : dans I'ensemble canonique, 
le nombre total d'impuretes presente des fluctuations d'ordre vN autour de sa valeur 
moyenne pN, alors que dans I'ensemble microcanonique, le nombre total d'impuretes 
est exactement pN dans tous les echantillons, et le desordre restant ne concerne que les 
positions des impuretes. II est clair que I'ensemble microcanonique est beaucoup moins 
desordonne que I'ensemble canonique, et en effet il conduit parfois a des reductions 
spectaculaires de bruit dans les etudes numeriques |147[[TK] . 

Pour la Chaine d'Ising avec couplages et champs transverses aleatoires, I'ensemble 
microcanonique au point critique ()11.2|) est defini par la contrainte globale |991 I66j 

L 
^(InJi-ln/ii) =0 (11.3) 

alors que dans I'ensemble canonique, la quantite du membre de gauche presente des fluc- 
tuations d'ordre \/L autour de sa valeur moyenne nulle. Les echantillons de I'ensemble 
microcanonique sont done 'plus proches' du point critique dans un certain sens, et 
permettent d'obtenir des resultats numeriques avec moins de fluctuations. 

Mais il y a un debat sur I'interet physique de cet ensemble microcanonique. D'un 
cote, certains auteurs |147M15j pensent que I'ensemble microcanonique doit etre prefere 
a I'ensemble canonique qui introduirait un bruit supplementaire susceptible de masquer 
les proprietes 'intrinseques' du systeme. D'un autre cote, les fluctuations d'ordre \/I 
d'une somme de I variables aleatoires sont precisement une propriete essentielle des 
systemes desordonnes, qui apparait par exemple dans les divers arguments physiques 
discutes dans I'lntroduction (PagelSJ. De plus, si on divise un systeme de taille L en deux 
sous-systemes de taille L/2, chaque moitie presente des fluctuations d'ordre \/L dans 
les deux ensembles : dans I'ensemble canonique, ces deux moities sont independantes, 
alors que dans I'ensemble microcanonique, les deux moities sont completement correlees 
et ont des fluctuations opposees. De ce point de vue, I'ensemble microcanonique pent 
done apparaitre comme artificiel. 

Bien sur, on s'attend a ce que les deux ensembles soient equivalents dans la limite 
thermodynamique, mais comme dans la theorie des lois d'echelle en taille finie, on 
mesure les exposants thermodynamiques sur les proprietes de taille finie, la discussion 
sur les deux ensembles conduit en fait au probleme de la theorie du 'finite-size scaling' 
pour les systemes desordonnes |147( I168( [T5] . 

11.1.3 Theorie du 'finite-size scaling' pour les systemes desordonnes ? 

En presence de desordre, la question centrale est la suivante |147l 11681 ITB] : si on 
obtient par simulations numeriques des valeurs pour une observable comme la sus- 
ceptibilite x{hL), pour plusieurs echantillons (i) pour diverses failles L, quelle est la 
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meilleure procedure pour analyser ces donnees ? La procedure usuelle consiste a moyen- 
ner I'observable sur les echantillons de taille fixee Xm{L) =< xih L,) >j et a analyser 
la dependance en L de la valeur moyenne Xm{L) comme dans la theorie du 'finite- 
size scaling' d'un systeme pur. Cependant, les references jl47( ll()8| IT^ qui ont etudie 
la question sont arrivees a la conclusion que la procedure correcte consiste a intro- 
duire une pseudo-temperature critique Tc{i,L) qui depend de Pechantillon. En effet, 
ces etudes ont montre que les fluctuations de xih ^) entre les echantillons venaient en 
grande partie des fluctuations de Tc{i,L), et que Temploi de la variable reduite de la 
variable (T — Tc{i,L)) dans I'analyse de 'finite-size scaling' permettait d'obtenir de bien 
meilleurs resultats (cf par exemple les Fig. 2 et Fig. 3 dans la reference |15j). De ce point 
de vue, I'avantage de I'ensemble microcanonique par rapport a I'ensemble canonique 
est justifie si la contrainte microcanonique determine en fait Tc{i, L) ou du moins reduit 
considerablement ses fluctuations entre echantillons, comme dans le cas des distribu- 
tions binaires J1471I15J . Pour les transitions de phase quantiques, la pseudo-temperature 
critique Tc{ii L) doit etre remplace par un pseudo point critique comme dans la chaine 
d'Ising avec couplages et champs transverses aleatoires, avec la contrainte (|11.3|1 . 

11.1.4 Observables interessantes en taille finie 

Les etudes numeriques sur les proprietes de taille finie dans les deux ensembles 
inni lEE] se sont interessees aux observables suivantes : 

• les aimantations de surface definies comme I'aimantation d'un bord lorsque le spin 
a I'autre extremite est fixe 

ml =< al > |^£=i (11.4) 

^< ai > 1,.=! (11.5) 



mj 



et leur correlation (?7ifr?T,£). 

• la correlation spin-spin entre les deux bords 

C{L) =< afal > (11.6) 

• le 'gap', qui represente la difference d' energies entre les deux premiers niveaux 

A(L) = Ei-Eo (11.7) 

Le but de ce chapitre est de comparer analytiquement les comportements de ces 
diverses observables dans les deux ensembles, au point critique et dans son voisinage. 

11.2 Resultats sur I'aimantation de surface 

11.2.1 Aimantation de surface et variable de Kesten 

L'aimantation de surface (|11.5|) dans un echantillon de taille (L + 1) a en fait une 
expression exacte en fonction des couplages aleatoires (HHl lES] 

mf = [1 + Zl]-i/' (11.8) 
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ou 

a une structure specifique de somme de produits de nombres aleatoires connue sous le 
nom de variable aleatoire de Kesten J114j . Ce type de variables apparait aussi dans les 
versions discretes du modele de Sinai J1581 llStil IH7] , et dans divers problemes formules 
en termes de produit de matrices de transfert aleatoires 2x2 |58l I35j . La version 
continue d'une variable de Kesten est la fonctionnelle exponentielle |22[ I136J 

rL 

Zl= / (ixe-^(^) (11.10) 

Jo 

du mouvement Brownien U{x) = Jq dyF{y) oii le processus {F(x)} correspond a la 
version continue des variables (— 21n(/ij/Jj)). 

11.2.2 Distributions asymptotiques de I'aimantation de surface 
11.2.2.1 Methode du col dans chaque echantillon 

II est naturel d'evaluer I'integrale (|11.1U|) pour L grand par une methode du col 

Zl ~ e^^ (11-11) 

L— »oo 

Oil {—El) < est le minimum du processus U{x) sur I'intervalle [0, L]. Au point critique, 
le scaling El ^ yL montre qu'il existe une distribution limite pour la variable d'echelle 
(In Zl)/VL ~ Inml/VL ~ El/VL. 

11.2.3 Lois d'echelle dans les deux ensembles 

Plus generalement, au voisinage du point critique on obtient les formes d'echelle 
suivantes pour la distribution de probabilite (— Inmf) 

2 ^/ -21nmf ^— \ 

PLi-lnml) = ^=Q{w = ^^;7 = ^V^ (11.12) 

VcrL V VcL J 

oil les fonctions d'echelle sont respectivement donnees dans les deux ensembles par 



Qcano{w,-f) =e{w>0) 



Vr J w-j ^TT 

2 



(11.13) 



Qmicro{w]j) =e{w>0)e{w>--f){2w + -f)e-'^'"-'" (11.14) 

Ces resultats out aussi ete obtenus dans la reference (^0] par une methode differente. 
L'avantage de la formulation par la methode du col (|11.11|) est d'une part de clarifier 
la validite des formules 1)11.14(1 et d'autre part, de pouvoir etre generalisee a d'autres 
observables, comme la correlation entre les deux aimantations de surface. 



Resultats sur Vaimantation de surface 93 

11.2.3.1 Distribution asymptotique de la correlation {mfrnf^) 

De meme, si on evalue par une methode du col dans chaque echantillon la correlation 
(mf m|^) entre les deux aimantations de surface, on obtient 

In(m^m^) ~ (11.15) 

Le calcul de la distribution asymptotique de (mfm£) se ramene done au calcul de la 
loi jointe de I'amplitude A{L) = Umax — Umin et du point final U{L) d'une trajectoire 
Brownienne partant de f/(0) = 0. Les resultats pour les deux ensembles (canonique et 
microcanonique) sont donnes dans la Publication [P15]. 

11.2.4 Moyenne de I'aimantation de surface 

11.2.4.1 Evenements rares 

On vient de voir qu'au point critique, I'aimantation de surface a pour comportement 
typique 

m{ = e-'^'" (11.16) 

oil w est une variable aleatoire d'ordre 1. La moyenne sur les echantillons de surface 
va done etre dominee par les echantillons rares qui presentent une aimantation anor- 
malement grande d'ordre 1, et sa decroissance avec la taille L va etre gouvernee par la 
dependance en L de la mesure de ces echantillons rares. Ainsi le comportement dans 
I'ensemble canonique [23 lEE] 



1 

I wv . 



l\]cano OC -j= (11.17) 



reflete directement la probabilite qu'une marche aleatoire ne repasse pas par I'origine 
dans I'intervalle [0, L]. En revanche, le comportement different dans I'ensemble micro- 
canonique |991 I66j 

[mWmicro OC - (11.18) 

pent etre interprete comme le rapport entre 

(i) la probabilite l/L^'"^ pour la probabilite de premier retour 
(ii) la probabilite 1/vX d'etre a I'origine en L. 

11.2.5 Resultats exacts par un calcul d'integrale de chemin 

Grace a une methode d'integrale de chemin avec contraintes (Publication [P15]), on 
pent calculer les distributions exactes de I'aimantation de surface en taille finie (c'est a 
dire sans faire la methode du col (|11.1H) ) et obtenir ainsi exactement les aimantations 
moyennes dans les deux ensembles. 
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Dans I'ensemble canonique, on trouve 

~ + 7 / dve~ 



[mlUnoil = I^V^) ~ ^'^^°^^"-' 



L^oo ^tigL 



+ 0(i) (11.19) 



alors que dans I'ensemble microcanonique, on trouve que les deux variables d'echelle 
sont differentes des deux cotes de la transition, avec les fonctions d'echelle suivantes 

K]™cro(7 = ^V^ > 0) ~ 1=^^ (11-20) 

L^oo \/aL 

K1.„.„„(P = ,^L < 0) ^.^ ^ ^'°° dpp-/^.-./,„ (2yf) K, (2y|.- j) + o (ij) 

Le regime critique du cote desordonne est done gouverne par la variable d'echelle 
p = [luL < dans I'ensemble microcanonique, alors qu'il est gouverne par la variable 
d'echelle 7 = [i\oL < dans I'ensemble canonique ! C'est parce qu'il existe en fait 
deux exposants de correlation. 

11.2.5.1 Discussion des deux exposants de correlation 

II apparait deux exposants de correlation distincts 1/ = 2 et P = 1 dans I'etude de la 
chaine dans la limite thermodynamique 75 : en particulier, la fonction de correlation 
moyenne est gouvernee par I'exposant z/ = 2, alors que la correlation typique est gou- 
vernee par I'exposant u = 1 [JS]- Plus generalement, ces deux echelles de longueur 
apparaissent aussi dans des modeles relies, en particulier dans I'etude des fonctions 
propres de I'operateur de Fokker-Planck du modele de Sinai (cf la discussion Page IHn|l 
et dans les modeles unidimensionnels de transport pour les fermions [51 ll39j . 

La presence de ces deux exposants pent etre explique comme suit [^. La premiere 



echelle de longueur correspond a la longueur ^ sur laquelle la valeur moyenne U{L) — U{{)) 
FqL = fiaL devient d'ordre 1, ce qui donne 

ir ^ avec i> = 1 (11.21) 

La seconde echelle de longueur correspond a la longueur S, sur laquelle la plupart des 
echantillons sont effectivement caracterises par une difference de potentiel {U{L) — U{0)) 
qui est du meme signe que la valeur moyenne : I'echelle vctL des fluctuations doit etre 
du meme ordre que la valeur moyenne, ce qui donne 

£ r avec u = 2 (11.22) 

Cette analyse montre qu'il n'est pas surprenant de trouver qu'une meme observable 
puisse etre gouvernee par deux exposants differents dans les deux ensembles : en effet, si 
on cherche la longueur sur laquelle la plupart des echantillons 'connaissent' vraiment le 
signe de Fq, on voit que c'est la longueur £ (|11.22|) dans I'ensemble canonique, mais que 
c'est la longueur £ (|11.21|) dans I'ensemble microcanonique. En conclusion, la definition 
meme des deux longueurs £ et £ montre que la contrainte microcanonique est susceptible 
de jouer un role dans les comportements d'echelle en taille finie. 
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11.3 Gap et correlation entre spins de bord dans les deux 
ensembles 

Contrairement aux aimantations de surface discutees jusqu'a present, le gap A(L) 
H11.7() et la correlation C{L) H11.6() entre les deux spins de bord ne peuvent pas etre 
exprimes simplement en termes de tous les couplages aleatoires de la chaine. Cependant, 
ils peuvent etre etudies j77j grace a la formulation en taille finie de la renormalisation de 
type Ma-Dasgupta, ce qui a permis de faire des predictions detaillees pour I'ensemble 
canonique [Hj, qui sont en excellent accord avec les donnees numeriques [ZZHnEl- Le but 
de cette section est de comparer ces resultats de I'ensemble canonique avec les resultats 
pour I'ensemble micro canonique calcules dans la Publication [P15]. 

11.3.1 Determination du gap et de la correlation par renormalisation 

L'analyse par renormalisation |751l77j permet de determiner le gap A(L) ()11.7|) et la 
correlation C{L) (|11.6j) dans chaque echantillon decrit par le potentiel Brownien U{x) 
for < X < L : le gap A(L) est asymptotiquement determine par la derniere barriere 
renormalisee ascendante G selon 

lnA(L) = -G (11.23) 

alors que la correlation C{L) est asymptotiquement donnee par la variable A = G — 
U{L) + U{Q) selon 

lnC7(L) = -A (11.24) 

En particulier, la renormalisation predit une relation tres simple dans chaque echantillon 
entre les deux observables 

G - A = [/(L) - C/(0) (11.25) 

Cette prediction a ete confirmee numeriquement au point critique, avec la convergence 
de la loi de la variable {G — A)/vL vers une distribution Gaussienne dans I'ensemble 
canonique, et vers une distribution 5 dans I'ensemble microcanonique _66j . 

11.3.2 Resultats au point critique 

11.3.2.1 Gap dans les deux ensembles 

Au point critique, les distributions de probabilite de la variable G = — InA(L) ont 
respectivement pour transformees de Laplace par rapport a la longueur L dans les deux 
ensembles 



+ 00 

dLe-P^ 



+ 00 

dLe-P^ 





sinh y/pG 
^ cosh ^/pG 



Pr\G) = ^— Vi--^^ (1126) 



p1^=°)(g) 



V^ ^ - VpG coth yTpG 



sinh y/pG 



g-v/p(jcotJi^p(j (11.27) 
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En particulier, les comportements en singularites essentielles a I'origine obtenus 
apres I'inversion des transforniees de Laplace, conduisent a des comportements en ex- 
ponentielles etirees pour les valeurs moyennes 

[AWir^r « LV6e-i(^^) ' (11.28) 

L—>oo 

[Aiirir = [Cl]™' « L-i/6e-i(2-'^)^/'+2(2-^i)^''« (11.29) 

11.3.2.2 Correlation dans les deux ensembles 

Dans I'ensemble canonique, la variable A = — InC(L) a une distribution tres simple 

P^aZ'oW = ^e-^ (11.30) 

qui conduit a une decroissance algebrique pour la valeur moyenne 

L—^oc Li 

alors que dans I'ensemble microcanonique, les variables G = A sont identiques : la 
distribution ()11.27|) est tres differente, et la valeur moyenne decroit selon I'exponentielle 
etiree (tTT^ . 

11.3.3 Comportements dans la region critique 

De meme, on pent comparer les resultats explicites sur le gap et la correlation au 
voisinage du point critique, du cote ordonne et du cote desordonne (Publication [P15]). 
On trouve que les distributions de probabilite sont gouvernees par I'exposant critique 
f = 2 dans les deux ensembles, mais que I'exposant u = \ gouverne les observables 
suivantes dans I'ensemble microcanonique, au lieu de I'exposant z^ = 2 de I'ensemble 
canonique : 

• la moyenne de la correlation dans la phase desordonnee 

• la moyenne du gap dans la phase ordonnee. 

De nouveau, ces differences entre les deux ensembles viennent de la presence de 
deux echelles de longueur (I11.21J) et (111.2211 . 

11.4 Conclusion 

Cette etude explicite des proprietes de taille finie de la transition de phase quan- 
tique, permet de clarifier I'origine des differences entre les deux ensembles, canonique 
et microcanonique, pour le desordre. Les variables d'echelle sont les memes dans les 
deux ensembles, mais les lois de probabilite reduites sont differentes, notamment dans 
leurs comportements asymptotiques. En consequence, les moyennes des observables qui 
sont gouvernees par des evenements rares peuvent avoir des comportements critiques 
tres differents. 
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Par ailleurs, I'etude de I'aimantation de surface qui admet une expression ferniee 
en ternies des variables de desordre, eclaire la signification de la methode de type Ma- 
Dasgupta. En effet, dans le cas particulier de I'aimantation de surface ecrite sous la 
forme de variable de Kesten, la methode de Ma-Dasgupta correspond exactement a 
la methode du col dans chaque echantillon. Plus generalement, on peut done voir la 
methode de Ma-Dasgupta pour les autres observables comme un acces direct 'au col' 
dans chaque echantillon, alors qu'on ne connait pas d'expressions fermees sur lesquelles 
on pourrait effectuer une methode du col usuelle. 

Publication associee [P15] 
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Conclusion 



L'interet des approches de type Ma-Dasgupta 

Dans ce memoire, nous avons decrit en detail comment differents types de systemes 
desordonnes pouvaient etre etudies par des procedures de renormalisation de type Ma- 
Dasgupta. L'interet de ces methodes vient des proprietes specifiques suivantes. 

• Leur domaine de validite : 

Elles permettent d'etudier les systemes desordonnes dans lesquels les heterogeneites 
spatiales du desordre dominent a grande echelle par rapport aux fluctuations ther- 
miques ou quantiques. Pour les systemes gouvernes par un point fixe de desordre infini, 
elles donnent des resultats exacts. Pour les systemes gouvernes par un point fixe de 
desordre fini suffisamment fort, elles constituent I'ordre dominant d'un developpement 
systematique par rapport a I'inverse du desordre. 

• Leur structure pour les regies de renormalisation : 

Contrairement aux methodes usuelles de renormalisation qui traitent I'espace de 
maniere homogene, les procedures de type Ma-Dasgupta sont inhomogenes dans I'es- 
pace pour mieux s'adapter aux realisations locales du desordre : c'est une variable 
extreme de desordre qui est decimee de maniere iterative et qui engendre un Hot de 
renormalisation pour une distribution de probabilite. En dimension d = 1, cette struc- 
ture particuliere permet d'obtenir des solutions explicites pour de nombreux modeles. 
En revanche, en dimension superieure d > 1, I'application des regies de renormalisa- 
tion n'a pu etre faite, du moins jusqu'a present, que de maniere numerique (cf Annexe). 

• Leur sens physique : 

Elles ont un sens physique tres clair, car les regies de renormalisation operent dans 
I'espace reel et sont basees sur des arguments physiques simples. Ainsi, dans les modeles 
classiques dynamiques, la renormalisation est basee sur le temps necessaire pour franchir 
une barriere de potentiel ou pour s'echapper d'un piege. Dans les modeles d'equilibre 
thermodynamique, comme les chaines de spins classiques ou I'heteropolymere a une 
interface, la renormalisation pent se voir comme un prolongement des arguments de 
type Imry-Ma sur I'energie aleatoire que le systeme est susceptible de gagner dans un 
domaine. Enfin, dans les modeles quantiques, la renormalisation permet de construire 
des amas de spins fortement correles. 
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• Leur resultats : 

Pour tous les modeles ou ces methodes peuvent etre appliquees, elles donnent une 
description tres complete des phenomenes. En effet, comme la procedure de renormal- 
isation est effectuee echantillon par ectiantillon, il est possible d'etudier par une meme 
analyse de tres nombreuses observables. Souvent, elles donnent meme acces aux dis- 
tributions de probabilite des observables par rapport a I'ensemble des echantillons, ce 
qui permet de bien comprendre I'influence des evenements rares sur certaines quantites 
physiques. 

• Leur efRcacite : 

Si elles fournissent des resultats qu'on ne pent pas obtenir par d'autres methodes, 
c'est finalement parce qu'elles reposent sur deux piliers tres differents, qui ne cohabitent 
generalement pas a I'interieur d'une meme approche, a savoir des arguments physiques 
heuristiques d'une part et des calculs exacts de renormalisation d'autre part. L'idee 
qualitative fructueuse des approches de type Ma-Dasgupta est d'integrer ces deux as- 
pects de maniere claire et coherente : (i) elles utilisent des arguments physiques simples 
pour identifier les degres de liberte importants a grande echelle, et pour definir une 
renormalisation directement sur ces degres de liberte (ii) elles etudient ensuite exacte- 
ment le Hot de cette renormalisation. 

Probabilites et renormalisation 

Pour conclure par une reflexion plus generale, la renormalisation est un langage 
naturel pour I'etude des probabilites. En effet, de meme que certains font de la prose 
sans le savoir, la theorie des probabilites a fait de la 'renormalisation' bien avant que 
les physiciens n'introduisent ce mot, avec le theoreme de la Limite Centrale, les notions 
de lois stables et de bassins d'attraction. 

Dans le cadre des systemes desordonnes, les questions probabilistes sous-jacentes 
relevent de la statistique des extremes [BH El] : I'equilibre a basse temperature est 
gouverne par la statistique des etats de basse energie ISHl IMl US OH \SM > alors que 
la dynamique a grand temps est gouvernee par la statistique des grandes barrieres 

[TniHTirTiiTTin]. 

De ce point de vue, les approches de type Ma-Dasgupta donnent un nouvel eclairage 
en proposant de maniere generale de renormaliser des distributions de probabilite par 
decimation locale d'une variable aleatoire extreme. En particulier, pour le cas des po- 
tentiels aleatoires unidimensionnels, elles proposent de calculer par renormalisation la 
distribution de probabilite jointe des positions des extrema et des barrieres qui les 
separent. Parmi les solutions explicites que nous avons decrites, la classe d'univer- 
salite la plus importante est la classe d'universalite des extrema Browniens, qui permet 
de comprendre de maniere unifiee I'equilibre ou la dynamique de differents types de 
systemes desordonnes. 



Annexe : Autres modeles etudies 
dans la litterature 



Comme nous I'avons deja mentionne dans I'introduction (PageEJ), les renormalisa- 
tions de type Ma-Dasgupta ont ete tres utilisees pour les systemes quantiques depuis 
les travaux de Daniel Fisher en 1994-1995 [7^1 EHI- Le but de cette Annexe est de men- 
tionner tres brievement les diverses directions qui ont ete suivies dans la litterature. 
(Des discussions generales sur la physique des systemes quantiques gouvernes par des 
points fixes de desordre infini se trouvent dans les references 



Diverses chames de spins quantiques 

La renormalisation de Ma-Dasgupta pour la chaine antiferromagnetique S = 1/2 
isotrope a ete generalisee dans les directions suivantes : 

• en presence d' anisotropic |Z01 

• en presence de dimerisation |97j 

• pour des spins plus eleves, en particulier 5 = 1 (cf ChapitreE3) et 5 = 3/2 jl49j . 



La reference j^S] contient une discussion du diagramme de phases des chaines antifer- 
romagnetiques en fonction des divers parametres anisotropie/spin/desordre. La renor- 
malisation a ete aussi utilisee pour caracteriser les proprietes de transport de ces chaines 

m- 

Parmi les autres types de chaines de spins etudiees, on pent citer 

• la chaine avec couplages et champ aleatoires jTSj (cf Chapitre [TT|) . 

• la chaine de Potts et la chaine 'clock' a q etats en champs aleatoires jl55j 

• le modele Ashkin- Teller ^6j 

• les echelles de spins jl32j . 

Mentionnons enfin pour terminer une etude recente sur un modele de bosons fS|. 
Cette liste non-exhaustive montre deja les nombreuses possibilites d'etude pour les 
modeles quantiques unidimensionnels. 
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Modeles quantiques en dimension superieure d > 1 

Pour les modeles quantiques antiferromagnetiques, la procedure de renormalisa- 
tion de Ma-Dasgupta a ete etudiee numeriquement depuis longtemps en dimension 
superieure ^Hl- L'etude recente [1 21] fait le point sur les points fixes obtenus en dimen- 
sions d = 2 et d = 3 pour le modele antiferromagnetique sur reseau hypercubique, ainsi 
que pour differentes variantes qui presentent de la dimerisation ou de la frustration. 

Pour le modele d'Ising avec couplages et champs magnetiques aleatoires (discute en 
dimension d = 1 dans le Chapitre lll|) . la renormalisation de type Ma-Dasgupta a ete 
etudie numeriquement en dimensions d = 2 et d = 3 J1381 175] . Le resultat essentiel est 
que la transition de phase quantique a temperature nulle est encore gouvernee par des 
points fixes de desordre infini en d = 2 et d = 3. 

Y a-t-il un espoir de calculer un jour exactement les exposants critiques en d = 2 ? 
D'un cote, la procedure de renormalisation de type Ma-Dasgupta detruit completement 
le reseau regulier initial qui devient un reseau aleatoire, et introduit de plus des correlations 
dans les variables de desordre, ce qui complique beaucoup I'analyse. D'un autre cote, la 
renormalisation de type Ma-Dasgupta transforme le probleme de la transition de phase 
quantique en un probleme de percolation de certains amas bidimensionnels |138j . Le 
calcul exact des exposants critiques en d = 2 sera done peut-etre possible un jour? 

Modeles classiques etudies dans la litterature 

Dans le cadre de la physique statistique classique, en dehors des modeles discutes 
dans le present memoire, les autres approches recentes de type Ma-Dasgupta concernent 
les problemes suivants : 

• la percolation 2D avec desordre constant par colonne |lU2j 

• le modele de Potts 2D dans la limite g — > oo |in3j 

• la transition de phase vers un etat absorbant en presence de desordre |lU4j 

• un modele d'exclusion en presence de desordre |in5j . 
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Resume 

Les procedures de renormalisation dans I'espace reel de type Ma-Dasgupta permet- 
tent d'etudier des systemes desordonnes gouvernes par des points fixes de fort desordre. 
Apres une presentation generale des idees physiques importantes et des methodes de 
calcul, ce memoire decrit les resultats explicites exacts que ces procedures de renor- 
malisation permettent d'obtenir dans difFerents niodeles unidimensionnels, classiques 
ou quantiques, dynamiques ou statiques. La majeure partie du memoire est consacree 
a des modeles de physique statistique, avec notamment 

(i) la dynamique hors equilibre d'une particule dans un potentiel Brownien ou dans un 
paysage de pieges aleatoires, 

(ii) la dynamique de croissance de domaines et I'equilibre thermodynamique des chaines 
de spins desordonnees classiques, 

(iii) la transition de delocalisation d'un polymere aleatoire a une interface. 
La derniere partie du memoire concerne deux chaines de spins quantiques desordonnees 
qui presentent une transition de phase a temperature nulle en fonction du desordre, a 
savoir 

(a) la chaine de spin S = 1 antiferromagnetique aleatoire 

(b) la chaine d'Ising avec couplages et champs transverses aleatoires. 

Abstract 

The Ma-Dasgupta real-space renormalization methods allow to study disordered 
systems which are governed by strong disorder fixed points. After a general introduction 
to the qualitative ideas and to the quantitative renormalization rules, we describe the 
explicit exact results that can be obtained in various one-dimensional models, either 
classical or quantum, either for dynamics or statics. The main part of this dissertation 
is devoted to statistical physics models, with special attention to 

(i) the off-equilibrium dynamics of a particle diffusing in a Brownian potential or in a 
trap landscape, 

(ii) the coarsening dynamics and the equilibrium of classical disordered spin chains, 
(iii) the delocalization transition of a random polymer at an interface. 
The last part of the dissertation deals with two disordered quantum spin chains which 
exhibit a zero-temperature phase transition as the disorder varies, namely 

(a) the random antiferromagnetic S = 1 spin chain, 

(b) the random transverse field Ising chain. 



